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Abstract 

This work is devoted to an examination of Stern-Gerlach forces consistent with special 
relativity and is motivated by recent interest in the relativistic Stern-Gerlach force acting 
on polarized protons in high-energy particle accelerators. The equations for the orbital and 
spin motion of a classical charged particle with arbitrary intrinsic magnetic dipole moment 
in an external electromagnetic field are considered and by imposing the constraints of 
special relativity and restricting to first order in spin (= first order h) a well-defined class 
of spin-orbit systems is obtained. All these systems can be treated on an equal footing 
including such prominent cases as those considered by Frenkel and by Good. 

The Frenkel case is considered in great detail because I show that this system is iden- 
tical with the one introduced by Derbenev and Kondratenko for studying spin motion 
in accelerators. In particular I prove that the spin-orbit system of Derbenev and Kon- 
dratenko is (nonmanifestly) Poincare covariant and identify the transformation properties 
of this system under the Poincare group. The Derbenev-Kondratenko Hamiltonian was 
originally proposed as a way to combine relativistic spin precession and the Lorentz force. 
The aforementioned findings now demonstrate that the Derbenev-Kondratenko Hamil- 
tonian also provides a legitimate framework for handling the relativistic Stern-Gerlach 
force. 

Numerical examples based on the Frenkel and Good cases for the HERA proton ring 
and electromagnetic traps are provided. 
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Introduction 



One of the most economical descriptions of the motion of polarized beams of spin 1/2 particles 
in high-energy accelerators is provided by the semiclassical Hamiltonian given by Derbenev and 
Kondratenko |pK73| | , which in the language of the Dirac equation can be written in the 
form 



Al,op AI ^opjOro Al, op, spin ' 



(0.1) 



with 



H 



AI ,op,orb 



/3 ■ i/vf^ -7?,, +m? + e-(h., , 



M^op^spzn op Ivl,op ' 



(0.2) 



where 



tTat = P«r — e ■ y4,, = canonical momentum vector 

A'l ,op J- M .op Ai .ov 

e 



W, 



"•"^ 2-m VJ, 
9-2 



[- \ 7^1 ■ P ■ ^At,op 



2 J,, ■ (J,, + 

AI,op \ AI,op ' 



2 ■ J,, J,, + m 

AI,op AI,op ' 



M,op M.op M,op 

^Af.op /\ ^M.op) ) + hermitian conjugate , (0.3) 



with 



J, 



M,op 



7rt 

AI.op 



TT., + 

AI.op ' 



(0.4) 



Furthermore I define 



M,op 



.if — if . ]-[ \ 

hi .op K'l .op AI.op AI.op/ ' 



^ AI.op — (-'^ AI.op ' '^M.op) 



-1/2 



I AI .op AI.op 



kinetic momentum vector 



mechanical momentum vector . 

(0.5) 



Here 



Aj^j Ej^j Bj^j represent the electromagnetic field, e, m denote charge and (non- 



vanishing) rest mass of the particle whereas g denotes the gyromagnetic factor [a]. The oper- 
ators I • a are Dirac matrices. This semiclassical Hamiltonian can be obtained from the 

' ' Ti op 

Dirac Hamiltonian (modified by the Pauli term) by a certain relativistic generalization of the 
Foldy-Wouthuysen transformation [6] in which terms of second and higher order in the spin are 
dropped. The orbital variables r^^ ^ ^Vm o ^^"^ spin variable a^^ obey the usual commuta- 
tion relations. The operator is often called the operator of 'canonical mean position', 

choose units where the vacuum hght velocity is c = 1. 
^Specifically r,^,^ tPmo obey the canonical commutation relations whereas 2 • t?^ /% obeys the same com- 
mutation relations as the Pauli matrices. 
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which explains the presence of the subscript 'M'. Note also that for vanishing electromagnetic 
field r,, _ is the position operator of Newton and Wigner [[NW49|| . The operator defines the 



M,op 



spin in a frame at rest w.r.t. the canonical mean position. The spin part H,, of the semi- 
classical Hamiltonian is in the form of a Stern-Gerlach (SG) energy term and the orbital part 
H,^ ^ resembles the standard orbital Hamiltonian. With this semiclassical Hamiltonian the 

M ,op,oro 

Heisenberg equations of motion for the orbital operators comprise the Lorentz force and SG like 
terms. The Heisenberg equation of motion for the spin operator is equivalent to the Thomas- 
Bargmann-Michel-Telegdi equation (Thomas-BMT equation) [pMT59| , |Tho27|| . This semiclas- 
sical Hamiltonian appears in the calculation of natural radiative polarization of electrons (the 
Sokolov-Ternov effect) ||DK73| , |Jac76| , ^T64|] in inhomogeneous fields. However, the interpre- 
tation of this Hamiltonian in terms of classical variables is also of great utility and has been 
used to describe the SG force for relativistic particles as well as to construct the action-angle 
variables for classical spin-orbit motion [ pHR94a| , [BHR94b| , PK73| , Per90a| , |Der90b| , |Yok87 



Since from the outset the semiclassical Hamiltonian is to be applied in a context where h can 
be neglected, the classical version contains the same information without additional approxi- 
mation. So it would be an unnecessary complication to continue to use the quantum version - 
one can just view the latter as a catalyst. Thus in the remainder of this paper I will use the 
classical form. 

My interest in relativistic SG forces stems from the suggestion that they could be used 
to separate spin states in storage rings and thereby provide polarized (anti-) proton beams 



951, ITTRETf 



Clearly, since the SG force is very small, one must be careful to include all 
terms in calculations. This is especially true in storage ring physics where there are not only 
strong transverse magnetic field gradients but also strong longitudinal gradients as well as large 
high frequency electric fields. One way to include all terms automatically is to use a formulation 
based on a Hamiltonian and of course one immediately thinks of the Derbenev-Kondratenko 
(DK) Hamiltonian. A formulation based on a Hamiltonian is also symplectic. The Hamiltonian 
approach has already been developed in ||BHK94a| , |BHK94b| ]. 

I have opened this paper with a review of the DK approach to spin and orbital motion but 
there is a large literature on classical, Poincare covariant equations of motion in the presence of 
SG forces going back as far as the paper of Frenkel of 1926 |[Fre26|| . See the reviews in pTSOj , 
|Nyb62| , |Pla66b| , [Roh72|, P?VW80| |. So, naturally, given the prominence of the DK Hamiltonian 
in the handling of spin effects in storage rings, one would like to know how it is related to other 
formalisms and whether it leads to Poincare covariant equations of motion. One would also 
like to know what general form the covariant equations of spin-orbit motion can take and to 
study the numerical implications of various choices. 

Although I also treat some other topics the main burden of this paper is then the following: 

• I prove that the equations of motion derived from the DK Hamiltonian are, after a 



transformation of the coordinates, identical to those parts of the Frenkel equations |[Fre2C: 
obtained after dropping second and higher order terms in the spin, and are therefore 
Poincare covariant. To the knowledge of the author the present proof is the first and it 
occupies most of the present work. Q The key to the proof is the realization that the 
position variable of the DK equations is not the spatial part of a space-time position 
but that this can be remedied if one first transforms it into the position variable given 
by Pryce which does not have this disadvantage. F| This leads to the conclusion that 



•^For the special case of the first order SG force a proof was effectively contained in [DS70|. 
■^I learnt this 'trick' from [DS70| who applied it to the special case of the first order SG force. 
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the approach of [ DK73 | is consistent with special relativity [c]. Thus the treatments in 
BHR94a| , PHR94b| , Per90a| , Per90b|| are also relativistic. 



• I then demonstrate, at first order in spin, how to construct more general classical Poincare 
covariant equations of spin-orbit motion and find that there is an unlimited number of 
possibilities allowed by the kinematic constraints with each possibility being characterized 
by the choice of five parameters. This unified approach enables one to classify equations 
appearing in the literature by identifying their characteristic parameters. In particular, 
besides the Frenkel equations, three other cases are considered in more detail: the Good- 
Nyborg-Rafanelli (=GNR) equations |Goo62, iNlyb64, liat7(J| | (which is also the approach 
chosen in |pPP95| , [NR87|| ), a case dealt with in ||GoM94|| , a case dealt with in ||Gos94|| , 
and the simple case of vanishing SG force. This technique could be extended to higher 
order in spin and the number of characteristic parameters would then increase. 

• Having achieved my main aims I then investigate in more detail another relationship 
between the DK equations and some of the other spin-orbit systems and show that some 
spin-orbit systems can be related by transformations. In particular I relate the Frenkel 
equations to the GNR equations. 

• I obtain numerical estimates of the SG force for simple spin and field configurations in 
the HERA proton storage ring ||Br95|| and in traps and compare the expectations from 
various equations. 

Even if one has a basis for choosing among the plethora of spin-orbit systems available, 
for example by experiment or by appealing to the Dirac equation, one must still decide which 
variables to use for representing the particle position. The natural choice is the Pryce position 
variable since it transforms like the spatial part of the space-time position. On the other 
hand, in Newton- Wigner coordinates the DK equations are in canonical form. As we shall see 
later, Newton- Wigner coordinates can lead to equations of motion very different from those 
for the Pryce coordinates. However, the Newton- Wigner coordinate differs by less than the 
Compton wave length from the Pryce coordinate, so that the relationship between the Frenkel 
and DK equations is not just mathematical but expresses the effective physical identity of the 
DK particle and the Frenkel particle. 

The paper is organized as follows.. In sections 1 and 2 the spin-orbit system defined by 
the DK equations is considered. In sections 3 and 4 I consider the spin-orbit system defined 
by the Frenkel equations and I show that it is the same system as in sections 1 and 2. The 
central part of the work is section 5 where the most general form for the spin-orbit systems 
under study is defined (equations (5.5)) and their main properties are explored. In sections 
6,7,8,9 I simply expose some further properties and give numerical examples. In particular in 
section 8 I consider the SG force in magnetic fields typical for HERA-p and in section 9 the SG 
force in electromagnetic traps. In section 10 I introduce the transformations between different 
spin-orbit systems, mentioned above. 

Throughout this work the time evolution of the spin is always that of | |BM'1'59| , |'l'ho27| | and 
the orbital motion is determined by the Lorentz force plus the SG force. The SG force consists 
of two parts, a part linear in the electromagnetic field vectors (the 'first order SG force') and a 
part quadratic in these fields (the 'second order SG force'). 
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1 The DK equations 
1.1 The DK Hamiltonian 

The classical DK Hamiltonian [d] is obtained from the semi classical version (0.1) by replacing 
the operators by classical canonical variables and the commutators by the following Poisson 
bracket relations ^ [e]: 



3 

Wj,(^k}M = J2 ^ikm- CTm ■ (j, = 1, 2, 3) (1.1) 



m=l 



Note that r^^ , t are the position and time variables used in the everyday business of acceler- 
ator physics if calculations are made in 'cartesian coordinates'. The electromagnetic field is 
characterized either by the potentials A^^^ or the field vectors E^^j, B^^. Thus one has [/]: 



dA 

B,, = V,, A A^^ , E,, = -V„0„ - . (1.2) 



The functions (p^j, Aj^_j^ E^^, B depend on r^^.t. This also fixes the meaning of the partial 
derivative d/dt in (1.2). The corresponding Maxwell equations are given in section 2. On 
introducing the quantities [g\: 



= PM-e-A 



M ' 



m 



2 



m VJ,,^ 2 J « 2 J,,-(^+m) - - - 
^ ^ ' '^M^Ej) , (1.3) 



2 ■ Jm Jm + 



the classical version of the DK Hamiltonian reads as [PK73|| : 



HA^m^Pm,^, t) = Jm + e ■ (t>M + ■ W,, . (1.4) 

To maintain consistency with the semiclassical nature of the DK Hamiltonian, throughout this 
paper all terms of second and higher order in the spin are neglected. This is the underlying 
approximation used in this paper. In particular one then always obtains the spin motion given 
by IIBMT59I, |rho27||. R 



1.2 The DK equations 

With this Hamiltonian the equations of motion for the orbital and spin variables - the DK 
equations - can be obtained using (1.1) and are: [] 

= R./.^aJa/ . (l-5a) 



^The variables f^^^'p^^,(j together with the Poisson bracket { , } define a Poisson algebra for the DK 
Hamiltonian. 

®This approximation is standard also in applications to accelerator physics. For the meaning of a power 
series expansion in spin, see also |Pla66b|. 

''in this paper the total time derivative is denoted by the prime '. 
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Pi, = {Pm, ^m}m = -r ■ (^M A BJ - e ■ V,,<P,, + — ■ (Tft^ ■ V,J/,, - V„((Tt . W,,) , 

SG term 

(1.5b) 

a'={a,i/„K, =#,,Aa, (1.5c) 



where the SG terms in the orbital equations have been exphcitly identified. The covariance of 
(1.5) under the Poincare group will be demonstrated in section 5. Note that (1.5c) is equivalent 
to the Thomas-BMT equation [ |BlVri'59| , [rho27|| . 
In the following it will be useful to define: 



M 



7 



9 -'2 7. 



■ ■ B • V 

7m + 1 " " " 



[- 



7. 



Af 



7m + 1 



]-(^M A^m) • 



By neglecting the spin terms in (1.5-6) one gets the familiar relations: 



(1.6a) 

:i.6b) 



;i.6c) 



^ ■ (7m ■ ^m)' = e ■ (v^, A BJ + e ■ , 
With (1.6) I find that (1.5c) reads as: 



7 



M ^ Af A-f 



CT 



;i-7) 



This is the usual form of the Thomas-BMT equation ||BlVri'59| , [l'ho27|| . The transformation 
properties of r„j,t, Vj^^, a under the Poincare group are discussed in detail in section 5. 



1.3 

By (1.1), (1.4), (1.5a) I have: Q 



fff ■ (^,, A E„) ■ 7f„ 



m ■ {Jj^j + m) 



■ (^M A ^) + 



e-g 



2 ■ m ■ J, 



■(^M Aa)- 



SG terms 

2 ■ m ■ J3 

Af 



SG terms 

0-2 /e 1 - R 

2 U J^,^ . ( J^,^ + m) - 



SG term 



e m + 2 ■ J 



!, ■ B,, ■ 7f + - 



TT 



m J3 . ( J^^ + m)2 « « m J,, ■ ( J^, + m) 



I ■ ^M ■ ^ • (1-8) 



SG terms 



^This equation is derived in Appendix B, see: (B.3-4). 
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Thus for the 'M' variables the canonical momentum is different from the kinetic momentum, 
i.e. 

which says that in the 'M' variables one has zitterbewegung. This effect disappears with the 
electromagnetic field. | 

The spin vector is constrained by: 

a^-a = h^/A. (1.9) 

Because (1.9) is of second order in spin it plays no role in this paper. It is only applied in sections 
8 and 9 for numerical calculations, where it gives the spin vector its correct length. Hence (1.9) 
acts as a numerical constraint to be inserted if the formulae are numerically evaluated. Note 
also that (1.9) is conserved under (1.5c). 

My next task will be to rewrite equations (1.5) in a form which facilitates comparison with 
other formalisms. 



Reexpressing the DK equations in terms of auxiliary 
variables 



2.1 



,Vj^,j,a by new variables rp,Vp,s which will later serve 



In this section I replace the variables f 
as the building blocks of the Poincare covariant formulae to be derived 
I define 



i 7 

f =f H ■ — — — -(aAv 



(2.1) 



P,op 



The quantity fp is the position variable of Pryce | )Fry49| ] . The corresponding operator 
describes (as t^^^^ does) a particle which is not pointlike from the point of view of the Dirac 
position operator. For the explicit form of these operators in the Dirac representation, see 
[Hei|| . For the special case of the first order SG force, see for example: [pS70| , PS72|| . Note 



that fj^j and differ by less than the Compton wave length. In contrast to f^^ , the variable fp 
transforms under Poincare transformations as the spatial part of the space-time position (see 
section 3). In the present paper, where I work classically, I do not rely on the Dirac equation 
but instead derive the Poincare covariance of the DK equations from the covariance property 
of fp. 

The components of fp are not canonical, i.e. 



{r 



r 

P,i 1 P,k 



3 

1=1 



— ■si + v 

7m 



M,i 



M 



7^0. (j. A; = 1,2, 3) (2.2) 



^The term 'zitterbewegung' was introduced in | ^ch3C | and in its original sense it only applies to a free 
quantum mechanical particle. Thus my term deals with the classical analogue and even applies in the presence 
of electromagnetic fields. For more details, see [ |Cor6S| , ^H62| 



^"The quantum mechanical proof of this covariance property, based on the Dirac equation (plus Pauli term), 
deals with the operator Vp^^ and is given for the special case of the first order SG force in | lDS72 |. For the 



general proof, see 
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But the complementary virtues of my two position variables are now clear; the variable is 
useful for studying covariance of the equations of motion (as seen below) whereas r^^ is useful 
for symplectic calculus (calculating spin-orbit transport maps) because the Poisson brackets 
for are canonical (see (l-l))- 

The velocity vector corresponding to fp reads as: 

V. ^ n - + ■ • A v^)) , (2.3a) 

and I define: 

7, ^ (1 - vl . Vp)-'/' 



(2.3b) 



Furthermore I define: 

s^7^-a--^-vl-a-v^, (2.4) 

from which follows: 



1 



. ■?! I . e . 



(7 = — • S + — - -V^-S-Vp , 

V ^ • (T = -U ^ • S , 

p p ' 

3 

{sj,Sk}M^J2^jki-{si + -fl-Vp^-vl-s) . (j, A; = 1,2, 3) (2.5) 

1=1 

One can therefore express the evolutions of T'^{t),Vj^{t), a{t) in terms of fp{t),Vp{t), s{t). The 
field equations for the 'M' fields are 

(2.6) 

which are the vacuum Maxwell equations. On introducing the abbreviations 

Ep = E^{fp,t) , Bp = B^{fp,t) , (2.7) 

one then obtains [h]: 

f)E f)H 
= VpABp--^ , = VpAEp + -^ , = Vl-Ep, = VI-Bp, 



(2.8) 



which also arc the vacuum Maxwell equations. The transformation properties of fp,t, Vp, s, Bp, Ep 
under the Poincare group are discussed in detail in section 5. 



2.2 

In Appendix A it is shown that the DK equations (1.5) lead by (2.4) to: 
s' = -. \^ -^p-vl-s- {Ep +v^ABp) + — -vl-Ep-s 



2/P p V — P ' -p- - — p/ 1 

7 



p 

p 



■^-^■l,-vl-s-v^^-Ep.Vp + ^-~.isABp)-^-.^.s^.Ep.Vp]. (2.9) 
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2.3 

Moreover in Appendix B it is shown that the DK equations (1.5) lead by (2.3) to: 

m • (7, • v^y = e • (i;, A 4) + e • 4 + ' ^ ■ V, (s^ -B^-El- {sA v^)) 
+|^-[2-s^-5;-t;,-5-(5At;jt.E;-^-, + (^-2)-(E;As) 



e2 



+ 77^ ■ [-(^ 2)' ■ < ■ 5p • (i?P A s) + {g- 2f -BI-B^- {v^ A s) 
-(^ - 2)2 . s t . 4 . 4 + (_^2 . . + 2 ■ (7 ■ t?; ■ t;^ + 4) ■ s ^ ■ 4 ■ 4 
-(^ - 2)2 • . 4 . (4 A s) + {g- 2)2 . (^7^ A B^Y ■ ■ {v^ A s) 
+{g - 2f . El ■ E^ ■ {v^ As) - {g' - i - g) ■ ■ B^ ■ {v, AB^) 
+{g-2).g-El-B^.s-(g-2)-g-v^^.E^-v^^-B^.s 
+(g-2)-g-v^^-E^-s^ ■B^-v^ + (g-2)-g-v^^-s-v^^-B^-E^ 
-{g'-i-g)-El.{v^A^-E^ + 2-g-El-{v^As)-{v,AB^)] . (2.10) 

2.4 

On collecting (2. 3a), (2. 9-10) one finds that the DK equations (1.5) are equivalent to: 

f;^v,, (2.11a) 
m • (7^ • v^y ^e-{v^AB^) + e-E^+ ' ^ ■ (s^ ■ B^ - eI ■ {sA v^)) 



V^-{2-s^-B'^-v,-g-{sAv,y-E'^-v, + {g-2)-{E'^As) 



e ■ 7 
"2^ 

+{g-2)-v^^-s-B'^-{g-2)-v^^-E'^-{v,As)] 

• [-(^ - 2)2 . t/t, . . {B^ A .?) + {g- 2)2 • fit . 5^ . (^^ ^ s) 

4r * 777- 

-(^ - 2)2 . s t ■E^-B^ + (-^2 . ^jt . + 2 . p . . + 4) • s t . fi^ . 4 
-(^-2)2.t;t .^.(4 As) + (^-2)2.(^^ A5,)t.i^.(^^;V^-^ 

+(^ - 2)2 -El-E,- {v, A s) - (/ - 4 • ^) • s t . B, ■ {v, A B,) 

+{g-2).g-El.B^-s-ig-2)-g.if'^.E^-if'^.B^-s 

+{g-2)-g-v^^-E,-s^ ■B,-v, + {g-2)-g-vl-s-v^^-B,-E, 

-(/ -A-g)-El-{v,As)-E^ + 2-g-El-{v,As)- {v, A B^)] , (2.11b) 

s'^^-[^-^,-vl-s-{E^+v,AB^) + ^-iP^-E^-s 



^ ^ lp-vl-^-^p-Ep-v, + ^--^-{sAB,)-^-.^.s^.E,.v,]. (2.11c) 

2 7p 27^ 



The covariance of (2.11) under the Poincare group will be demonstrated in section 5. By 
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neglecting the spin terms in (2.11b) one gets the famihar relations: 

m-{j^-vj = e-{v^AB^)+e-E^ , 7; = ^ . t/t . . (2.12) 



The equivalence between the DK equations and the 
Frenkel equations 



3.1 



In this section the DK equations are shown to be equivalent to the Frenkel equations. 

P P 

To study the covariance properties of , {Tp , s*, 1 1 first construct the variables X , U where 
the components are defined by 



< = (r1,^■^) 



(/i = l,...,4) 



(3.1) 



and represent the space-time position. The components denote the corresponding 4- velocity, 
i.e. □ 

K =K^Yr = ■ 1 < = ■ ^p' ^ ■ ^-)^ ' = 1' 4) (3.2) 

where r denotes the proper time. 

The utility of the spin vector s becomes apparent by encoding it in a spin tensor S . To 
achieve this I first introduce the 'dipole moment tensor' 



M. 



IIP 



/ 







^^p,2 


i 


^p.i 









f^p.i 


i 


^ P,2 











i 


^P.S 


V 


—I ■ 6p -^ 


—I ■ Ep^ 










{fi,u = l,...,4) (3.3) 



which by definition transforms as a tensor of rank 2. One calls flp the 'magnetic dipole moment' 
and Ep the 'electric dipole moment' [ |[\lyb64| |. The magnetic dipole moment resp. the electric 
dipole moment in the rest frame is denoted by (!„ resp. and by definition it is the intrinsic 

p 

magnetic dipole moment resp. intrinsic electric dipole moment of the particle. The tensor M 
appears frequently in the literature on the relativistic SG force and it occurs already in the 
1926 paper of Frenkel |[Fre26|] . Note that M is also used in the theory of relativistic fluids 
resp. composite particles. Q 

Although M is especially useful for treating particles with combined intrinsic magnetic 
and intrinsic electric dipole moments, such particles are not treated in this paper. In this study 
I only consider particles without intrinsic electric dipole moment. Denoting the dipole moment 
tensor in the rest frame by M , I then have: 



M. 



/ 







^^R,2 











^^R,l 

















v 











/ 



(/i,z/= 1,...,4) 



^^In this paper the upper dot symbol denotes the total proper time derivative. 
^^See for example the textbook treatments in |DS72, Moe72, Syn58|. 
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Bearing in mind that for a charged particle the intrinsic magnetic dipole moment is related to 
the rest frame spin by: 



^ _ e-g ^ 



(3.4) 



one obtains: 



M, 





1 





0-3 


-(72 


M 


e-g 




-0-3 





0"! 





2 • m 




0-2 


-CTi 










[ 











; 



(/i,z/ = 1,...,4) 



j-> 

One now defines the spin tensor 5* as that tensor, which in the rest frame reads as: 

/ ^3 -^2 \ 



s 



-(T3 












-(Ti 








\ 








J 



(/i,i/ = 1,..,4) 



(3.5) 



Thus the spin tensor S is given by: 



M = ■ S . 



2 ■ m 



By transforming back from the rest frame (see subsection 7.2) one obtains: 

(/i,j/ = l,...,4) (3.6) 



S, 



( ° 


S3 


-S2 


-i-qi \ 


-S3 





Sl 


-i ■ q2 




-Sl 





-i ■ q3 




i ■ q2 


i ■ q3 


J 



q = s AVj: 



so that one has 



e-g 
fip = s , 



e-g 

2 • m 



■ i^P A s) 



(3.7) 



^ 2-m 

For a particle with vanishing intrinsic electric dipole moment equations (3. 2), (3. 6) enshrine the 
kinematic constraints [k]: 



p p 

u, -u^ =-1, 

<.-f/"=0. (/i = l,...,4) 



(3.8a) 
(3.8b) 



With (3.6) I have encoded the spin vector s in the spin tensor S . For the consistency of the 

p 

rank 2 tensor property of S with the transformation properties of the spin vector s and of Vp, 
see subsections 5.4 and 5.5. 

I also introduce the rank two antisymmetric tensor field describing the electromagnetic field 
and whose components are defined by 








Bp,. 


-^P,2 


—i 


■Ep., 


\ 









^P,l 


—i 


■ Ep^2 






Bp,2 


-^p,i 





—i 


■ Ep^3 




[ 


^ ■ Ep, 


i - 2 









/ 
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The corresponding operator S is obtained in [FG61a, HW63]. 



(/i,Z/ 



.4) 



(3.9) 
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I.e. 



Then [/] with the definition: 



d 



d d (9 , 



p /tJ. 5 



the vacuum Maxwell equations (2.8) read as: 



= 



(/i = l,...,4) 



(3.10) 



(/i,z/,p = l,...,4) (3.11) 



3.2 

I can now reexpress the evolution (2.11) of rp(t),Vp(t),s(t) in terms of X^ (t) , (t) , (r) . 
In subsection 5.2 it is shown that equations (2.11) are equivalent to: 

C/; , (3.12a) 



m ^ 



4 ■ 



P P P P P P P P 



+ 



(9 



. '-^-s^^-k-u^-o^f: 

2 ■ ^ 



e^-i9-'2)-g p p p p 
e-g 



2-m 
2 ■ m 



p p p p 



p p p p 

Spuj ■ F^x -Ux -U^ 



4 ■ m'^ 
4 • 



p p p p 

SuLO ■ ^wX -Ux 



' pv 



u: 



UJp 



u: 



p p 

^Xw ' ^LUX ) 



(3.12b) 



,4) (3.12c) 



These are the equations given by Frenkel in 1926 |Pre26|| . |^ They respect the kinematic 
constraints (3.8) and I can therefore conclude that: 

• The DK equations are equivalent to the Frenkel equations. 

In this derivation it was essential that the electromagnetic field obeys the vacuum Maxwell 
equations. The covariance of (3.12) under the Poincare group will be demonstrated in section 



5. Note also that (3.12c) is equivalent to the BMT equation ||BMT59| , [Nyb64|| . The term 
e-g 



2 ■ m 



p p 
F ■ S* 



P P 



p p p p p p p p 

Sfuv ■ Fivx -Ux -U^ + S^^ ■ F^x -Ux -U^ 



is that part of the rhs of (3.12c) which is independent of the forces acting on the orbital motion. 
The remaining term 



p p 



u: ■ u: 



represents the Thomas precession [ l'ho27 ]. 

It follows from the normalization of the spin vector that: 

^pi^ ' ^fiu = ^ /2 . 

As in (1.9) this equation is of second order in spin so that it plays no role in this paper. Note 
also that (3.13) is conserved under (3.12c). 



(3.13) 



^■^Equations (3.12) are equivalent to equations (13a), (14), (21), (21a-b) in |Fre26|. They were rederived by 
many authors. See also the reviews in | |BT80| , |Nyb62| , |Pla66b| , |Roh72| , [rVW80| . 
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4 Rederiving the Frenkel equations in terms of a Hamil- 
tonian which is a Poincare scalar 



4.1 

The reader will perhaps be interested to learn that one can rederive the Frenkel equations from 
a Hamiltonian with the proper time as independent variable. 

P P P 

One begins by noting that X , P and S obey the following Poisson bracket relations 



P P 

P P P P P P 

{x;,x:}, = {x;,5:,}, = {p;,p;}, = {p;,s:j, = o. {^^,u,\u = l,...A) (4.1) 

Then I introduce the abbreviations 



4 ■ m 



M =m+ S^^ ■ F^^ , 



n; ^ p; - e ■ a; = (v?;, < )^ , = 1, 4) (4.2) 

where Ap{rp,t), (f)p{rp,t) are the potentials of the electromagnetic field so that: 

^ ^ ^ ^ ^ BA 

Bp = VpAAp, Ep = -Vp(t)p - -rrr ■ (4.3) 



i.e: 



dt 



<. = «-^<- (/i,z/ = l,...,4) (4.4) 



The Hamiltonian is [m]: 

^" = (--^)-n:-n: + i-M" + ^ii^-n"-5" -Pl-n:, (4.5) 

and the corresponding equations of motion are: 

= {Xfi }p , 

P r^^P ^^P 



U, ={U,,H }p = {{X^,H }p,H h 



^. = {<.,^"}p- (/i,z/ = l,...,4) (4.6) 



By evaluating the Poisson brackets one finds that these are the Frenkel equations (3.12). 

Because S is a tensor of rank 2 and F is a tensor field of rank 2 it follows by (4.2) that 

P P P 

M is a scalar field. Moreover, because P is a 4- vector and A a 4- vector field, one finds that 
all three terms of the Hamiltonian (4.5) are scalar fields. 

^^In the quantum mechanical analogue to the system (4.1) the position operator correponding to must 
be different from fp because this system involves (unlike the system for the 'M' variables) negative energy 
states. See for example [JM63, Cor68|. 
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4.2 

With (4.5) one gets: 



m ^ 4 ■ m ^ I - 

4 ■ m'^ 



+ P^,-K^-K-Fl-Sl^. (/i = l,...,4) (4.7) 



and: 



With (4.8) one can write (4.7) in the elegant form: 



n; = e ■ f;, . t/; - ^ . ■ , (/. = i, 4) (4.9) 



which will be useful in subsection 7.3. 

Using (4.8) the constraints (3.8), expressed in terms of the variables X ,P and S , read 

as: 

n'' ■ n"^ = - 2 ■ m ■ m"" , 
5;,-n" = o, (/i = i,...,4) 



so that the Hamiltonian if vanishes, if the constraints are taken into account. One thus has 
a constrained Hamiltonian system. The constraints are to be taken into account only in the 
final results (e.g. the equations of motion). 

4.3 

In the 'P' variables the canonical momentum is vf^ and using (4. 2), (4. 8) and Appendix D I get: 

m-'yp-Vp = 7fp^ r ■ ■ 5^ ■ TT^ - Ei- {sA vr^) ■ vf^ 

m"^ 2 ■ m'^ ■ 7p 

e-(g — 2) ^ e ■ (q — 2) , -± . 

+ ^ 2 ■ • ^ ■ + ■ 7p ■ {E^ A s) . 

1 ■ m'^ 2 ■ m 

1 thus observe for the 'P' variables that the canonical momentum is different from the kinetic 
momentum, i.e. 

As in case of the 'M' variables (see the end of section 1), one has zitterbewegung for the 'P' 
variables. This effects disappears with the electromagnetic field. Q 

^^One can modify the Hamiltonian (4.5) in a way such that zitterbewegung arises even for the free particle. 
Then higher orders of the spin are important iCor6^ , Fia66b| . 
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5 Relating the Frenkel equations to other approaches. 
The Poincare covariance 



5.1 

The Frenkel equations are just a special case allowed by the kinematic constraints (3.8) and 
in fact various forms for the SG forces are possible even if one requires, as I do, that the 
spin equation is equivalent to the BMT equation and that the electromagnetic field obeys the 
vacuum Maxwell equations. |^ In fact the constraints (3.8) allow the following generalization 
of (3.12): 

K = K ' (5-la) 
^ = ^-<.-^r+<, (5.1b) 
• p ^ e ■ g ( p p _ p p 

2 Tfl I '^'^ u)v '-^fj.uj uju 



e-{9-2) / P 



2 ■ m 



■,p ^^p ^^p ^p ^p ^^p ^^p 



S,.-F^,-U, -U, -S,^-F,,.U, -U , (/x,z/ = l,...,4) (5.1c) 



P P 

where the 4-vector Y collects the SG force terms. I require that Y only depends on the 

P P P P 

following dimensional quantities: m, U , e ■ F , S , d . If one assumes that the dependence 
on U ,e ■ F , S ,d is polynomial and in particular of first order in 5* , then by dimensional 

P P 

analysis the dependence on F resp. d is at most quadratic resp. linear and the constraints 
(3.8) lead to the following most general ansatz [n]: 



^P e ■ C2 ( r^P r^P r^P , TT^ tt^ t-^^ 

4 ■ m' 



~^ ^ o ^2 ^ ■ '^io ■ Ux ■ dxF^^ + ^ ■ is^^ ■ d^F + U ■ S^^ ■ ■ d^F^^ 



I ^-^7 _ ^P ,r.Pj^P ■ C3 . . . jf _ • ^4 _ ^P _ jjP T^P 

2 ■ m? ^ A - rrfi A-m? ^ 

+ ^ ■ ■ ■ ^p"- ■ ' = ^' (^-2) 

where ci,...,C7 are dimensionless real numbers. This can be further simplified by using the 
vacuum Maxwell equations (3.11), so that 

Sl ■ <C = -2 • C ■ . (/i = 1, 4) (5.3) 

Hence by (3.11) the terms on the rhs of (5.2) which are proportional to Cg, Cj are not independent 
from the others so that the general ansatz (5.2) finally simplifies to: 

= ~ 4 . ■ {^^^^ ' ^P^uju + Uf, ■ S^^ -Ux ■ dx F^^ 

e ■ (C2 - ci - 2) P P P P P e'-C3 . . . . 
2 ■ rri'^ ^ A ■ m-^ ^ p ^ 

■ ■ u: ■ s:, ■ i^pi + ^ ■ ■ ■ ■ f/: • (/^ = i, 4) (5.4) 



,p ^p „p ^^p 
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An early work about the nommiqueness of the relativistic SG force is: [ Moe49[ . 
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Combining this with (5.1) I get the following general system of equations: 



X 



U. 



^ m ^ 



e ■ c-i 
4 ■ 



p p p p p p p 



e • (c2 - ci 



2-m? 



5, 



e 



2-m 
e-{9-2) 
2 ■ m 



p p 
F ■ 9 



p p 

Q ■ F 

flUJ LOU 



P P P P 



P P P P 

Syu ■ - Ux - 



(5.5a) 



(5.5b) 



(/x,i/ = 1,...,4) (5.5c) 



The covariance of (5.5) under the Poincare group will be demonstrated in subsection 5.3. 
The Frenkel equations correspond to 



Cl 







C2 = g , C3 = {g-2) ■ g , = g 



C5 



ig-2)' 



(5.6) 



and by neglecting the second order SG terms in the Frenkel equations one gets the 'reduced' 
Frenkel equations, which correspond to: 



Cl = , C2 = g , C3 = C4 = C5 = . 
The GNR equations, which are defined by [poo62| , [Nyb64| , [Raf70|| , correspond to: 

Cl = g - 2 , C2 = g , cs = C4 = C5 = . 
Another interesting choice is given by | |Cos94|| : 

Cl = -2 , C2 = g , C3 = g^ , 04 = g , c^ = -g"^ + 2- g 



This is a modification of a force given in [|CoM94|| by a 'redshift term': 



Cl = -2 
A very simple choice is 



C2 = g 



Cl 



cz = g 



C2 = C3 



C4 = 



C5 



C4 



C5 =0 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



corresponding to a SG force, which vanishes for the 'P' variables. 

That the kinematic constraints (3.8) allow the general form (5.4) of Y illustrates that the 
five parameters ci, ...,C5 are, on the classical level, just phenomenological constants which can 
only be fixed by comparison with experiments (e.g. in storage rings with polarized beams) [o]. 
The physical implications of this plethora of possibilities will be addressed later (see sections 
8,9). Note that in the above derivation of Y the assumption of first order dependence on spin 
was essential. The reader who is interested in higher order spin terms is advised to consult the 
large literature on classical relativistic spin. See, for example, the book ||Uor68|j for references. 



A recent interesting treatment, nonlinear in spin, can be found in ||Cos94 



"'^^Note also that for every choice of the characteristic parameters Ci,...,C5 the following equations of the 
previous sections remain valid: (1.1-4), (1.5c), (1.6-7), (1.9), (2.1-9), (2.11a), (2.11c), (3.1-11), (3.13). 
"'^^See equation [4.57] therein. Note that I neglect second order spin terms. 
^"See equation [15] in [CoM94|. Note that I neglect second order spin terms. 
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5.2 

It follows from (5.5b) and by using Appendix D that: 

m • (7^ • Vp)' ^ e - {Vp ABp) + e- Ep + ^ ^ " • (s^ ■ - El ■ (s A 

2 ■ m ■ V 

+^ • [(ci + 2) ■ . 5; ■ t;, - C2 • {sA v^y ■ 4 • Vp + (c2 - Ci - 2) ■ {E'^ A 5) 

+(C2 -c,-2)-v^^-s-B'^- (C2 - ci - 2) • • • (^7, A s)] 

+ IT-i ■h-vl-Bp-{BpAs)-c,-Bl-Bp-{vpAs) + C5-s''-Ep-Bp 

+ (-C3 ■ t;^ ■ - C5 + 2 ■ C4 - C3) ■ ■ 5^ ■ + C5 ■ ^/^ ■ • (ip A s) 

-C5 ■ {Vp A Bp Y ■ Ep ■ {Vp A s) - C5 ■ # ■ ■ {Vp A s) 

+ (2 ■ C4 - C3) ■ ■ 5p ■ (^Tp A Bp) + c^-El-Bp- s- ■ v^^ ■ E^ ■ v^^ ■ Bp ■ s 

+C3 -v^^-Ep-s^ -Bp-Vp + c^-v^p-s-v^^-Bp- Ep 

+(2 • C4 - C3) • • (iTp A • 4 + 2 • C4 • • (iTp A • (iTp A Bp)] . 

Collecting this with (2.3a), (2.9) I get 

r'^^Vp, (5.12a) 

m • (7p • Vp)' = e-{VpABp) + e-Ep + j^;^^ ■ Vp (s^ ■ Bp - El ■ {sA Vp)^ 

■ [(ci + 2) ■ . . - C2 ■ (s A VpY ■E'^-Vp + (C2 - ci - 2) • (E; a s) 

2 ■ m 



+(C2 - ci - 2) • • s • 5; - (C2 - ci - 2) • • • {vp A s)] 

+ ■[c,-v^^-Bp-{BpA^-c,-Bl-Bp-{vpAs) + C5-s''-Ep-Bp 

+ (-C3 ■ wj, ■ Wp - C5 + 2 ■ C4 - C3) ■ ■ Bp ■ Ep + C5 ■ ;7], • Ep ■ {Ep A s) 
-C5 ■ {Vp A Bp)^ ■ Ep ■ {Vp A s) -c^-El-Ep- {Vp A s) 
+ (2 ■ C4 - C3) ■ ■ Bp ■ {Vp A Bp) + C3 -El -Bp- s- ■ ■ Ep ■ vl ■ Bp ■ s 
+C3 ■ ■ Ep ■ s'^ ■ Bp ■ Vp + ■ vl ■ s ■ ■ Bp ■ Ep 

+(2 ■Ci-c^)-El-{VpAs)-Ep + 2-Ci-El-{VpAs)- {Vp A Bp)] , (5.12b) 
s'^^-[^-jp-vl-s.{Ep + VpABp) + ^.v^^.Ep.s 

III Z 

-3—^.^^.vl-s-iP^-Ep.Vp + ^-.^.{sABp)-^-.^.s^.Ep.Vp]. (5.12c) 

z ^ Ip ^ Ip 

It is shown in Appendix D that (5.12) is equivalent to (5.5). For the special choice (5.6) this 
means that (2.11) is equivalent to (3.12), i.e. that the DK equations are equivalent to the 
Frenkel equations! 

5.3 

In this subsection I show that equations (5.5) are Poincare covariant. Note that transforms 

P P P 

as a space-time position, U as a 4-vector, S as a tensor of rank 2 and F as a tensor field 
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of rank 2. Thus (5.5) are covariant under the restricted Poincare group [p] because equations 
(5.5) are constructed from covariant quantities. To show that (5.5) is covariant under the 
whole Poincare group I first consider space inversion (=parity transformation). This acts in 
the following way: 



r — »• r 



d d d d 



dXl ' dX^ ' 9X3" ' dX, 



4 



Bp ^ Bp , 
Ep -Ep , 
s — > s , 

q^-q. (/i=l,...,4) (5.13) 



Hence one sees that with the building blocks: 



no pseudotensors occur in (5.5), so that (5.5) is covariant under the parity transformation. 
With (5.13) one sees that rp,Vp, Ep,q are polar vectors whereas Bp,s aie axial vectors. 

Next I consider the time inversion (=time reversal transformation). This acts in the follow- 
ing way: 



r ^ -r , 

x;^ifi,-x;),, 



p 

tJ, 'V I P " p ' " I p J 1 

d d d d 



t/r (-7p 



^ "-dxrox^dxr dxl'^' 

Bp — > -Bp , 
Ep^ Ep , 
s ^ —s, 

q^q. (^ = 1,...,4) (5.14) 

Now I use the fact that (5.12) is equivalent to (5.5). In fact, the application of (5.14) to (5.12b- 
c) shows first of all that the spatial parts of (5.5a-c) are covariant under the time reversal 
transformation. It then follows by using (3.8) to get 

f/f = - ■ ( f/f ■ u; + ■ + u; -u^), 

q = sAVp^ sA{VpABp + Ep) vl ■ Ep ■ q , (5.15) 

m m ' 

that the temporal parts of (5.5a-c) are covariant under the time reversal transformation. Hence 
one concludes that (5.5) is covariant under the time reversal transformation. In summary: 
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Note that the direct product or the contraction of two tensors is again a tensor, not a pseudotensor. For 



the distinction between tensors and pseudotensors, see for example [Moe72 
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• Equations (5.5) are covariant under the Poincare group. Specifically, they are invari- 
ant under Poincare transformations (except that the electromagnetic field transforms as 
described in section 3). 

• As a special case the Frenkel equations (3.12) are Poincare covariant from which it follows 
that also (2.11) is Poincare covariant [q]. 

Because transforms as a tensor field of rank 2 one also finds that (3.11) is invariant under 
the Poincare group. 



5.4 

P 

So far I have assumed that S is a tensor. Now for completeness and as announced in section 3 

p 

I demonstrate that the tensor property of S is consistent with the transformation properties 
of Vp under the Poincare group, i.e. the relation q = s A Vp is conserved under Poincare 
transformations. In fact, by (5.13-14) one sees that it is conserved under spatial rotations, 
space inversion and time reversal and in the remaining part of this subsection I show that it is 
also conserved under the proper Lorentz transformation. 

The infinitesimal proper Lorentz transformation (^infinitesimal Lorentz boost) is defined 

by 

J. boost ~ 
J. boost 

^ boost _ boost 

L:r = l, {j,k = 1,2,3) (5.16) 
where the infinitesimal vector v, ^ denotes the relative velocity of the frames connected by 

boost •J •J 

- boost 



L . From (5.16) one obtains 

boost boost _ , ^ J \ /- «\ 

V -^^P -^-P^ (i/,p = l,...,4) (5.17) 

which is consistent with the fact that ij""""* belongs to the Lorentz group. The spin tensor 
transforms under L via: 

_P „P,boost boost boost _,P , ^ j\ / — 

Sau^S^, =L„. -L,^ -S^. (/x,i/ = 1,...,4) (5.18) 



Abbreviating: 



^Pjboost / _^ . \ 

^ ^ V^boost' ^ ^boost ) ' 



one sees by (5.18) that s, q transform under iJ'"""'' via: 

s^s. = s — iJ, , A (f , 

boost boost ' 

^ ^ ?,oost ^loost ^ S . (5.19) 

Because C/^ is a 4- vector, one observes that Vp,^p transform under i^*""""* via: 

v„ ^ v„, . = v„ + v[ ■ v„ ■ v„ — V, , . (5.20) 

P Pjboost P boost P P boost \ / 
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Collecting (5.19-20) one observes: 



boost boost Jr*. boost ' 

so that the relation: q = sA Vp is conserved under the infinitesimal proper Lorentz transforma- 
tion L . This concludes the proof that this relation is conserved under the whole Poincare 
group. 

5.5 

Subsection 5.3 listed the transformation properties of rp,t,Vp, s, Bp, under the Poincare 
group. Combining these with (2.1), (2. 3a), (2. 4-5), (2. 7) one also obtains the transformation 
properties of r^jyt, Vj^^, a, Bj^^, Ej^^. One observes: 

• The relations (2.1), (2. 3a), (2. 4-5), (2. 7) are invariant under Poincare transformations. 

• By (2.7) the 'P' fields depend on rp,t in the same way as the 'M' fields depend on Vj^j^t. 

Because (5.12) is Poincare covariant, one concludes that also the DK equations (1.5) are 
Poincare covariant. Specifically, they are invariant under Poincare transformations (except 
that the electromagnetic 'M' fields transform in the same way as the 'P' fields) [r]. 

Note that ifl^^i ■ t)^ does not transform as a space-time position and (7J^^ • v^^^i ■ 7^)^* is 
not a 4- vector but both transform nonlinearly in a complicated way. Analogously one obtains 
the well known property that the rest frame spin vector a is not the spatial part of a 4- vector. 
Using (2.5) and the tensor property of I now discuss the transformation properties of a. 
Firstly by (5.13), a is, like s, an axial vector and by (5.14) it transforms in the same way under 
time reversal. 

Secondly, under the infinitesimal proper Lorentz transformation (5.16) it transforms via: 

^^^+^^-K.o.At;,)Aa, (5.21) 
7p + 1 

which is simply a rotation associated with a change of orientation of the reference frame due 
to the boost. Thus a has in fact the transformation properties of the rest frame spin vector 



BMT59, Jac75, Tho27 . 



5.6 

As already mentioned, r^j, t are the position and time variables used in the everyday business 
of accelerator physics. In fact via (2.1) and by neglecting spin: = fp, so that in this 
approximation 

(flft), (/i = l,...,4) (5.22) 
transforms as a space-time position. Also one observes in this approximation that: 

(Im-vI^-Im), (/i = l,-,4) (5.23) 

is a 4-vector. Even if one includes the spin one obtains these transformation properties if the 
SG-force is neglected. In fact in that approximation one gets: 

^;=^M, (5.24a) 
^ ■ (7m ■ ^mY = e ■ {v,, A BJ +e-E^^, (5.24b) 
a' = n^,^Aa . (5.24c) 

21 



BMT59 
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6 Using the spin pseudo-4- vector 
6.1 

So far I have described the spin in terms of the spin tensor S but because the particle only 

has a intrinsic magnetic dipole moment and no intrinsic electric dipole moment one can also 
describe the spin just in terms of a pseudo-4-vector T defined by: 

^/T = ■ ^M^P- ■ sip -Ul , (/i = 1, 4) (6.1) 

where e^upi^ is the Levi-Civita symbol. From this it follows that 

= -i ■ e^upL. ■ Tp -Ul . (/i, = 1, 4) (6.2) 
The constraints (3.8) now read as: 

K ■ K = -1 ' 

t; ■ f/; = . (6.3) 

From (2.4-5), (3.2), (3.6) and (6.1) one has: 

Tp={flX),, (/i = l,-,4) 

fp = — -s + 'yp-vl-s-Vp = a^ ■ vl ■ a ■ , 

T4 = i ■ ■ fp = i ■ 'jp ■ vl ■ s = i ■ ■ vl ■ a . (6.4) 

Thus under the space inversion one gets by (5. 13), (6.1): 

T;^(ft,-T;)^, (/i = l,...,4) (6.5) 

so that in fact is a pseudo-4- vector. Under the time reversal one gets by (5. 14), (6.1): 

T;^(-f;,T;),. (/x = i,...,4) (6.6) 

6.2 

Now I introduce the pseudotensor field F dual to F defined by 

Fpu = ^ ■ ^^lyp^ ■ . (/^, = 1, 4) (6.7) 



For more details on t"", see for example |BMT59| , |Cor6S| , |FG614 |Nyb64|] 



Note that e is the totally antisymmetric pseudotensor of rank 4 with £1234 = 1. 
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from which it follows that: 



In Appendix E it is shown that (5.5) is equivalent to: 



K = K ' (6-8a) 
j-jP _ e p p e-C2 p p p ~p e ■ (ci + 2) p p p p p ~ p 



2 ■ ^ 4 ■ m'^ ^ m 

2 2 

6 ' i"* ^ p p 6 ' ^ p p p p p 

5- ■ Pau ■ Fya " " 1 ? " -^ui' '^u '^a '"^0 ' ^aB ) (6.8b) 



^ ' ^ e ■ (^ - 2; _p „p „p „p 

2 ■ m ■ ■ " 2 ■ m 



^ = f-!: ■ • < + • ■ T: . Fl . U: . (/. = !, 4) (6.8c) 



Note that (6.8) conserves the kinematic constraints (6.3) and (6.8c) is the BMT equation 
BMT59|| . Note also that (6.8) is covariant under the Poincare group [^because (5.5) is, too. 



It follows from the normalization of the spin vector that: 

K ■ K = n'/A . (6.9) 



As in (1.9) this equation is of second order in spin so that it plays no role in this paper. Note 
also that (6.9) is conserved under (6.8c). 



7 The nonrelativistic limit. The rest frame 
7.1 

In this first subsection I consider the nonrelativistic limit (= zeroth order in 1/c) and in the 
remainder of this section I consider the particle rest frame. I do this for the general case, i.e. 
for arbitrary values of ci, ...,05. In the nonrelativistic limit (5.12) leads to: 

f'=Vp, (7.1a) 



p ^p ' 
m ■ V 



^ 2 ■ m 2 ■ m at 



22 2 

+ ^ .El-B^-s + . ( -C3 + 2 ■ C4 - C5 ) ■ ?t . 4 . + ^ ■ s t -E^-B^, 

4 ■ 4 ■ 4 ■ 

(7.1b) 



One now sees that the second order SG terms survive even in the nonrelativistic limit. This 
should come as no surprise: second order terms are seen in the usual semi-relativistic Foldy- 
Wouthuysen transformations [b]. Note that the third term on the rhs of (7.1b) is sometimes 

^H.e. (6.8) is invariant under Poincare transformations (except that the electromagnetic field transforms in 
the prescribed way) 

^^The partial derivative d/dt in (7.1), (7.2) and (7.5) acts on functions depending on rp,t. 
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called the 'magnetodynamic force'. In the nonrelativistic limit the Frenkel equations lead via 
(5.6), (7.1) to: 



,„.;T;=e.4 + f^-V,(.rt.B,)-l^.(?A^) 
^ 2 ■ m 2 ■ m ot 

4 ■ m2 p p ^2 p p 4-^2 p p ' 

.-" = ^-(sA5,), (7.2) 
2 ■ m 

which is equivalent to: |^ 

----e.E,, + ^.V,.(at.i^J-^.(.A^) 
eM^ + 2)^ R - e'-^?' z? R 

^' = |:^-(^A5J, (7.3) 
In the nonrelativistic limit the GNR equations lead via (5. 8), (7.1) to: 



r 

p 



m ■ V 



The choice (5.9) leads in the nonrelativistic limit via (7.1) to: 

^ e ■ q -i, . e ■ q ,^ dE„ , 

^ 2 ■ m 2 ■ m ot 

4 ■ ^ 4 ■ 

s' = ^-{sAB,), (7.5) 
2 ■ m 



which agrees with [4.57] in [|Cos94|| . For the choice (5.11) one gets in the nonrelativistic limit 
via (7.1): 

m ■ = e ■ Ep , 

To summarize the above one observes that the characteristic parameters C3, — C3 + 2-04 — 05, C5 
only vanish if = C3 = C4 = C5. Using simple linear algebra it follows from (7.1b) that 

^^The partial derivative d/dt in (7.3) acts on functions depending on , t. Note also that in the nonrelativistic 
limit one has: s = a. 
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knowledge of the nonrelativistic limit (7.1) uniquely determines the values of the characteristic 



parameters Ci 
(see eq 



C5. As an immediate application I note that the particle described in ||CPP95 
1.1) thereof) has a nonrelativistic limit identical to (7.4). If this particle belongs to 
the class obeying (5.5), one can conclude that the particle described in ||CPP95|| obeys the GNR 
equations, i.e. the characteristic parameters ci, ...,05 for this particle assume the values given 
by (5.8). Thus ||CPP95|| provides a good example of how to identify Ci, C5. That in fact the 
particle described in | |UPP95| is characterized by the choice (5.8) of the parameters Ci, ...,05 is 
supported by the results of section 8. 



7.2 

Although the behaviour in the nonrelativistic limit is closely related to that in the rest frame, 
which I define as the frame for which the velocity Vp vanishes, the two behaviours are neverthe- 
less distinct. In the nonrelativistic limit I am observing the motion in a selected inertial frame 
but motion with respect to the rest frame is actually motion with respect to an accelerated 
frame. The Lorentz transformation L to the particle rest frame transforms X to the rest 
frame space-time position X^: 



X. 



(/i 



.4) 



(7.7) 



where 



R 
R 

4i 

R 
AA 



(7p - 1) 



If 



Z ■ 'Jp ■ V 

-L 



P ^P,3 ' 
R 

jA ! 



V^-Vp 



(j,fc = l,2,3) 



6. 



jk 



Since the rest frame is an accelerated frame, L depends on the proper time. Relative to the 
rest frame the particle motion vanishes, i.e. 



U. 



(0,^) 



(^=1,...,4) 



(7.9) 



In the rest frame the quantities S , F transform to: 



jP qP 



/ 







— 0"2 




-0-3 





0"! 











V 














/ 



...,<) 



f7.10a) 



(/i,i/=l,...,4) (7.10b) 



where I also used (3.5). Note that on combining (7. 8), (7. 10a) one obtains (3.6). Furthermore 
one concludes from (6. 4), (7. 8): 

T; = L;,.T; = (crt,0)^, (/i = l,...,4) (7.11) 

where denotes the spin pseudo-4-vector w.r.t. the rest frame. 
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If N'^ (r) transforms under the Lorentz group as a 4- vector or as a space-time position and 



N^^{t) as a tensor of rank 2, the rest frame proper time derivative {d/dr)^ is defined by 



,^).<. = ^^A ■ C ■ |: (^A.-^ ■ L^,-' ■ K,) . {^,,u = 1, 4) (7.12) 



The distinction between {d/dT)„ and d/dr takes into account the proper time dependence of 

p 

the rest frame which occurs because in general the 4-velocity U is not a constant of motion. 
Applying (7.12) to the quantities X ,U ,S yields: 



.J^)rSIu = lIp-LI^-SI^. (/i,i^ = 1,...,4) (7.13) 



Combining (5.1) with (7.13) results in 



)nX, = f/, , (7.14a) 



, d ^ ^^R ^^R 

=^-c- < ■ < + ■ , (7.14b) 

( d . ^R _ ^R jR ( e-g p p _ p p 

' ■ ~ [SL ■ Kx ■ Ul ■ U; - S: ■ . Ul ■ U:]] , {f,,u = l, 4) (7.14c) 



2 ■ m 



where in (7.14a) I also used (7.9). To simplify (7.14) I introduce the abbreviation 

W' = ^^-^^^ 

from which it follows by (7.10b) that: |^ 

^IKp = ■ ■ ■ dlF^p . (/i, z/, p = 1, 4) (7.16) 
Combining (5. 4), (7.10), (7. 16) one gets 

e ■ (C2 - Cl - 2) R R R R R ■ C3 R R R R 

■ ■ ■ + ' ' ^'^ ■ ' P 

■ Ku-U^ ■ ■ F + Ji ■ S ■ F^ ■ F -U^ , (/i = 1, 4) (7.17) 



4-m3 " "^P P^ A-m? 



''The partial derivatives always act on functions depending on Xi , X2 , ^3 , 
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so that (7.14) specifies to: 

. d 



U., 



e ■ (c2 - ci - ^ 



^ 4 ■ m2 



u, 



sl ■ f/f 



■ C4 



2 ■ 

R 



R R 



■ U- 



n« p« , e2 ■ C3 



4 ■ m'^ ^ 



4 ■ m3 



pLU ^ UJ 1 



'dr'"' 
e-{9-2) 
2 ■ m 



e- g 

2 ■ m 



pw 



5, 



R 



s 



R 



pa; 



ojp 



s 



pa; 



^i-i^^A-f/r-f/; 



(7.18a) 



(7.18b) 



. (/i,z/ = l,...,4) (7.18c) 



Note that (7.18) has the same form as (5.5). Thus the use of (d/dr)^ provides an economic for- 
mulation of the rest frame equations of motion. To compare the rest frame equations (7.18) with 
the nonrelativistic behaviour studied in the previous subsection I introduce the abbreviations 

K = (^l ^4 )m ' ^p = (7« ■ ?^!,, ^ ■ 7fl)p , 

f"^(4,-2-4). (/i = l,...,4) (7.19) 



The simple structures of U and S allow the spatial parts of equations (7.18a-c) to be easily 
obtained and one gets: ^ 

i-^)R^^R = , (7.20a) 



, d , ^ ^ e ■ Co ^ ^ . i ■ e ■ (co — 2 — ci) ,^ dE 

m- (^)„v„ = e-E„ + - VJo^ ■ B) ^— ^ • a A 

^dr'"" « ^ 2-m 2-m ^ dXf 



R 



2 2 2 

+ ^—^ ■Ei-B^-a + ■ (-C3 + 2 ■ C4 - C5) ■ • 5^ ■ + -a^ -E^-B^, 

4 . ^2 i-m^ A-rn? r r ^ 

(7.20b) 
(7.20c) 

By (7.20a) one sees that the rest frame is that frame where the derivative {d/dr)^ of the position 
vanishes. Note that (7.20) also displays the fact that the rest frame is an accelerated frame. 
In fact, by the definition of {d/dr)^ one sees that the rest frame proper time derivative of a 
vanishing quantity in general is nonvanishing, which is exemplified by (7.20a-b): 

7.3 

In the special case of the Frenkel equations it is interesting to consider the rest frame behaviour 

p p 

of n instead of U . Abbreviating: 

n; ^ lI, ■ n: ^ (vf^, n,"), , (^ = i, 4) (7.21) 
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The nabla operator V^^ acts on functions depending on , ~i ■ , a and it is the gradient w.r.t. r^. 
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one obtains by (4. 9), (7.12): 



, d ■, R 6 R R G ' CI ^R ^R R / . \ 

— )„n„ =-.F„„-U„ -^^■S„..,-d„F. (/i = l,...,4) 



m 



4 ■ ^ 



Together with (7.20a), (7.20c) one then gets: 

(|:)»^-» = 0, 

d 



(^)A = -4 + |^-v.(.-*.B„). 

, (i , ^ e ■ q , ^ , 
(—)„(?= ■ (aAB^) . 



(7.22) 

(7.23a) 
(7.23b) 
(7.23c) 



Comparing with (7.20) one sees that in the rest frame the canonical momentum vector for the 
Frenkel equations fulfills the same equation of motion as the kinetic momentum vector of the 



GNR equations. Note that (7.23) was obtained also in [pia66a|| 



8 Estimating the strength of the SG force in magnets 



8.1 



This paper is partly motivated by the suggestion ||CPP95|, [NR87|| that the SG force can be used 



to separate spin states, either in real space or 'energy space', in (anti-)proton storage rings. 
Now that I have general forms for the relativistic SG force I am in a position to carry this 
study further. In this section I will apply my equations of motion to the HERA proton ring 
(HERA-p) | |Br95| . I do this for the general case, i.e. for arbitrary values of Ci, C5. Note that 



for (anti-)protons one has: {g — 2)/2 ^ 1.79. 



8.2 

In this subsection I study the equations of motion of the Pryce coordinates. I will only 
the case of static (i.e. time independent) magnetic fields and vanishing electric fields, i.e. 



~dt 







. 



I leave it to the reader to investigate other cases. To facilitate the numerical evaluations 
use Gaussian units |^ so that I drop the convention: c = 1. In the case of a static maj 
field (5.12b) reduces to: 



treat 
: □ 

(8.1) 

I will 

cnetic 



(C + 



e ■ C2 



2 ■ m ■ c ■ 7^ 



2 ■ m ■ 

^2 



■ [(ci + 2) ■ ■B'^-v, + (C2 - ci - 2) ■ ■ s- 5; 
■ [C5 -vl-B^- {B^ As}-C5- Bi- B^- {v^ A s) 



~'~ 4 ■ ■ 

+ (2 ■ C4 - C3) -s^ -Bp - {Vp A B 



(8.2) 



^^The partial derivative d/dt in (8.1) acts on functions depending on f^jt. 
■^"See for example [Jac75|. 



28 



To simplify the comparison with | CPP95 ], I express sin terms of a in (8.2). This leads by (2.4) 
to: 31 

6 -» 



e ■ C2 



2 ■ m ■ c ■ 7^ 



7; 



2 ■ m ■ 



[(ci + 2)-(7,-a- 



c2 ■ (7p + 1) 
c2.(7, + l) ""-^ ^ 



■ ^2 . (^"^^ 1) ■ ^ ^ ■ ^p) -C5-lp- Bl- Bp - {Vp A <j) 



+ (2-C4-C3)-# ■(7p-a 



7 



■ (7p + 1) 



-Vp ■ Vp) ■ {Vp A Bp)] 



- ■ (VpABp) + 



e ■ C2 



2 ■ m ■ c ■ 7p 



■Wphp-Bl- 



a 



7; 



+^-[(c. + 2).(7...^ 
+ {c2-c,-2)-v^-a-B' 



7 



c2 ■ (7p + 1) 



— ■ CT^ ■ V ■ B^ ■ V 

c'-ilp + l) ^ - - - 
■ ■ Vp ■ Vp)"^ ■ B'^ ■ Vp 



7 



C2 ■ (7p + 1) 



, . % — 7 ■ \c5 -ip-v^ - Bp - (Bp A 

[-C5 + 2 ■ C4, - C3) ■ a'' ■ Vp ■ vi ■ Bp ■ {Bp A Vp) 



-C5 --fp-Bl-Bp- {Vp A <t) + (2 ■ C4 - C3) ■ 7p • 5t . a • {Vp A Bp)] 



- ■ C A5J + - ^ X . 



2 ■ m ■ c 2 ■ m ■ 



-(ci + 2) 



+7p ■ (C2 



7, 



7^ 



rt 



„ , X v ■ a ■ v'^ ■ B' ■ v„ 
c2-(7^ + l) ^ P p p 



7p 



7p + l 



c^-2)-vi-a-B'] + 



[^l-{c, + 2)-a^ -B'^-Vp 



[C5 --ip-vl-Bp- {Bp A a) 



4 ■ m? ■ 



[-C5 + 2-Ci-C3)-vl-a-vl-Bp- {Bp A Vp) 



,2 . (^^ + 1) V - ■ - - -^P - P p 

-C5 -^p- B^-Bp- {Vp A a) + (2 ■ C4 - C3) --/p-a^ -Bp- {Vp A fi^)] 
I say that the magnetic field is 'transverse' if 



(8.3) 



^p-Bp = 
B'=0. 



(8.4) 



Hence for a static, transverse magnetic field (8.3) simplifies to: 

^ ■ (7p ■ Vp)' = - ■ {Vp A Bp) + ■ #,(at . 5 ) 

c 2 ■ m ■ c 



31r 



The nabla operator Vj, in this paper always acts on functions depending on rp , i, , <? and it is the gradient 
w.r.t. Tp. 

■^^This corresponds, for example, to the case of a particle travelling instantaneously parallel to the axis in a 
pure quadrupole magnetic field or to a particle travelling instantaneously perpendicular to a pure dipole field. 
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4 ■ ■ 



■ [-C5 ■ 7p • 5t . fi^ . (t;^ A a) + (2 ■ C4 - C3) ■ 7p ■ ■ ■ (t^p A B^)] . (8.5) 



If the magnetic field is transverse and the spin is parallel to the magnetic field then (8.5) 
simplifies to: 

^ ■ (7p ■ v,y =--{v,A 4) + ■ ■ 4) 

c 2 ■ m ■ c 

(,2 ^ ^ 

+ 2. 4 ■ ■ + 2 ■ C4 - C3) • ■ 5^ • (y^ A B^) . (8.6) 

Thus with C2 = g and when using Pryce coordinates the first order SG force equals the SG 
force given in ||CPP95|] . Note that the SG force in (8.6) does not depend on Ci. 
For protons one has 

= 10-23 . ergs ■ Gauss'^ = 6.3 • 10"^^ ■ MeV ■ Gauss-^ . (8.7) 

m ■ c 

In the HERA proton ring at about 800 GeV the magnetic field Bhd in the bending magnets 
is about 45 kGauss and the quadrupole magnetic field gradient Ghq is about 9 kGauss/cm. 
Thus one has 



e^-Blr,-h _ 

cm 



- 6.5 ■ 10-2^ ■ dyne = 4.1 ■ 10"^^ ■ MeV —'^ 



m? ■ 



e ■ (^HQ ■ ^ _ 9 1 . 10-20 . ^y^^ ^ 5 7 . 10-14 . MeV ■ cnT^ . (8.8) 



m ■ c 



For the static, transverse magnetic field with spin vector a parallel to the magnetic field the 
numerical values (8.8) lead via (8.6) to the following maximal values of the SG force: 

^'^^ •Vp(a^-Bp) ^ C2 ■ 2.3 ■ 10"2° ■ c/yne = C2 ■ 1.4 ■ 10"^^ ■ MeV^ ■ cm-i 



2 ■ m ■ C 

-C5 + 2 ■ C4 - C3) --ip-a^ - Bp - {Vp A B^ 



(-C5 + 2 ■ C4 - C3) ■ 7^ ■ 8.2 ■ IQ-^^ ■ dyne = (-C5 + 2 ■ C4 - C3) ■ 7^ ■ 5.1 ■ 10"^^ ■ MeV ■ cm'^ 
= (-C5 + 2 ■ C4 - C3) ■ 7.3 • 10^^^ ■ dyne = (-C5 + 2 ■ C4 - C3) ■ 4.6 ■ 10"^^ ■ MeV ■ cm'^ , 

(8.9) 

where I assumed 7p = 900. 

Thus at HERA-p energies and for the Frenkel equations (where one has: —05 + 2-04 — 03 = 4) 
(see (5.6)), the second order SG force in the dipoles completely outweighs the first order SG 
force in the quadrupoles. It is also simple to show that the second order force for a particle 
travelling about 1 mm off axis through a quadrupole is negligible compared to both of the 
above. In the case of the GNR force (see (5.8)) the second order force is zero. Clearly, if the SG 
force is to be utilized in high energy proton storage rings, one must first decide which equation 
of motion appertains. 

I say that the magnetic field is 'longitudinal' if 

Vp ABp = , 

B'p = ■ vl -B'^-Vp, (8.10) 

%-Vp 
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so that for a static, longitudinal magnetic field one gets: 



• v'^ ■ B • V 

p 



Hence for a static, longitudinal magnetic field (8.3) simplifies to: 

e ■ C2 1 



m ■ (7^ ■ Vp)' 



2 ■ m ■ c ■ V, 



■ ■ a ■ B' . 



(8.11) 



Note that the second order SG terms have disappeared and that the force only depends on C2. 
If the spin is parallel to the field this also agrees with the form given in | |CPP95| | if C2 = g. 
However, because Ci does not appear in (8.6) and (8.11), one sees that for these special magnetic 
fields, the GNR and reduced Frenkel equations are identical. 

I consider a gradient Giong of the longitudinal field along the longitudinal direction of about 



9 kGauss/cm ||CPP95|] , so that: 

6 ■ Grlong ' ^ 



m ■ C 



9.1 ■ 10"^° ■ dyne = 5.7 ■ 10"^^ ■ MeV ■ cm'^ 



(8.12) 



Thus for the static, longitudinal magnetic field with spin vector o parallel to the magnetic field 
the numerical value (8.12) leads via (8.11) to the following maximal value of the SG force: 



e • C2 



2 • m ■ c v^^ ■ Vp 



^p-vl-3-B'^ 



C2-lp- 2.3 ■ 10^2° ■ dyne 



C2 ■ 7p ■ 1.4 ■ 10"^^ ■ MeV ■ cm'^ = C2 ■ 2 ■ 10"^^ ■ dyne = C2 ■ 1.3 ■ 10"" ■ MeV ■ cm'^ 



(8.13) 



This agrees with [ UPP95 | where it is pointed out that the longitudinal SG force is much larger 
than the first order transverse force. It is also much larger than the second order transverse 
force. However, before this force could be used to separate the spin ensemble into two parts 
|CPP95|| , a way must be found to overcome the severe mixing ||Hof96|| caused by incoherent 



synchrotron oscillations. For a survey of other problems see |Per95| , Per90b |. 

Because for the GNR equations C2 = g, the identity of (8. 6), (8. 11) with the corresponding 
results in ||GPP95|| is consistent with the supposition in section 7 that the particle described 
in ||CPP95|| obeys the GNR equations, i.e. belongs to the choice (5.8) of the characteristic 
parameters Ci, C5. 



8.3 



This paper opened with a description of the DK Hamiltonian. This is based on the 'M' variables 
and so far I have only used these to provide the missing link connecting the DK equations 
to special relativity. However, since many other investigations |PHR94a| , PHR94b] , Per904 



Uer90b | have been based on this Hamiltonian it is natural that one inspects the equations of 
motion for the 'M' variables in more detail. As in the previous subsection I will only treat the 
case of static (i.e. time independent) magnetic fields and vanishing electric fields, i.e.: Q 



dt 











(8.14) 
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The partial derivative d/dt in (8.14) acts on functions depending on r,^j,t. 
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Note that by (2.7) the conditions (8.1), (8. 14) are equivalent. 
I begin by using Appendix C to rewrite (8.3) as: 

^ • (7m • = -■{v^^ BJ + . (C2 - 2 + -^) • V^(at ■ BJ 

c z ■ rn ■ c 7J^^ 

+^-r4T-[-^-(^-2)-ci-7M]-^^--^M-t^:-^M--^M 

+(o - 2) • -v^ -B' -B] 



e 

4 



7m 

2 



4 • • 



([4-A.^^+^^.{2.c,-Cs)]-a^-B^-{v^ABj 



+4 • • [-2 • C4 + C3 + C5 + 4 + (^=^ - 2 • y) • -^—r] ■ vl ■ a ■ vj, ■ B.. • (v,^ A B 



r2 ^ _L 1 L ^ ' -J ' I - I ^ , -|J -M ^ K^M'^^MJ 

+ [(^-2)-^--^-7^-C5]-4-S^-(aAS^)-C5-7M-^l-^M-(^M Aa) 

2 

+(^-2)^-^^-(^Ai;,)t.5^.5^ 

7m + 1 

-;^-(/-2-^)--^-(5'Ai;^)t.S^.i;i.S^.O ■ (8-15) 

If the magnetic field is transverse and the spin is parallel to the magnetic field then (8.15) 
simplifies to: 

^ • (7m ■ ^m)' = ; ■ (^^m a B^) + - . (C2 - 2 + -^) • V^(a^ • B^) 

+ , 4 -[4-4-7M+7M-(2-C4-C3-C5)]-at.^^.(i;^A^J , (8.16) 

which for the Frenkel equations (see (5.6)) becomes: 

^ • (7m • <y = 7 • (^M A BJ + . (^ - 2 + . v^(at . BJ 

+;;^ ■ • ^M • (^^M A BJ . (8.17) 
For the GNR case (see (5.8)) equation (8.16) becomes: 

^ • (7m • O' = ; • (^M A BJ + . (^ - 2 + -^) • V^(at . 

c z ■ m ■ c 7j^^ 

• (7m - 1) • • ^m • (^M A S^) . (8.18) 



The first order parts of (8.17) and (8.18) are identical. 

Although I have obtained the equations for the kinetic momentum, Hamiltonians lead more 
naturally to equations of motion for canonical momenta. If the magnetic field is transverse and 
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the spin is parallel to the magnetic field then by (B.4),(B.6),(B.17),(B.25-26) and by taking 
just the first order SG terms the Frenkel case gives: 

2-m-c 7^, 



{9-2) 



l^.v^ 

2 ■ m ■ c3 7^,, + 1 



a-B'. 



(8.19) 



and 



Pa 



m ■ c 



A 5. 



■{9-2 + 



+ 



2 , 

7m' 



V, 



5, 



7. 



e ■ - 2) 

2 ■ m ■ 7^^ + 1 



■ o 



B' 



(8.20) 



In equations (8.15-20) the first order SG piece contains the factor {9/2 — 1 + 1/7^^-) which differs 
from the corresponding term for 'P' variables in (8.6) by the term (—1 + 1/7^^) [|Hof95 1 . This 
term can be directly traced back to a similar term in Vtj^^ in (1.6c) associated with Thomas 
precession [|Jac76|] . "^^ 

Further differences appear for the second order terms both for the Frenkel and GNR equa- 
tions. Thus naive use of the DK Hamiltonian to obtain estimates of the relativistic SG force 
could lead one to the conclusion that it gives a first order SG force very different from that 
predicted by ||CPP95|| . In particular, for electrons one has: {9 — 2)/2 0.00116 so that the 
Thomas term causes a massive relative change. 

However, it is clear that the difference is only an artifact of the choice of position coordinates 
|PH95|| : the rest frame implied by the 'M' variables is different from that of the 'P' variables 
and the corresponding Thomas precession terms are different. Furthermore, although the forces 
on Tj^j and Vp can be rather different, the position variable r^^ is always closer to Tp than the 
particle Compton wave length (see (2.1)). Of course, the corresponding equations of motion 
conserve this property in time. 

The equation of motion (2.36) given in [ !JJKP96 | also contains the (—1 + 1/7^^) term and 
thereby appears to differ from ||CPP95|| . This should now come as no surprise since the former 
works with Newton- Wigner coordinates. 



9 Estimating the strength of the SG force in electromag- 
netic traps 

9.1 

Most accounts of the SG force emphasize the first order component. But now that one has 
seen that the second order force can be important at high energy in large fields it is interesting 
to estimate their effect under other circumstances where the SG forces play a central role. An 
obvious case is that of the nonrelativistic SG force in electromagnetic traps [pS V86|| . As in the 
previous section I use Gaussian units and I consider the general case, i.e. arbitrary values of 
Ci,...,C5. By (7.1b) one has: 

^ e ■ C2 ^ ^ . e ■ (co — 2 — ci) ,^ dE„. 
m-vL = e- Ep + - ^—.S/p{a^ -Bp) \^ ^ . ^ ^^ 



p 



2-m-c ^ ^ 2 - m - dt 



•^^Note that in the nonrelativistic hmit there is no difference (see (7.2-3) or put jj^^ — 1) 
■^^The partial derivative d/dt in (9.1) and (9.2) acts on functions depending on rp,i. 
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If the magnetic field and electric field are parallel to the spin then (9.1) simplifies to: 

2-ci) ,^ dE, 



^ 2 ■ m ■ c 2 ■ m ■ at 



■ El-B.-a. (9.2) 



2 ■ ■ ^ 

9.2 

To apply my equations to the case of the electron I calculate 
e ■ h 



m ■ c 



1.9 ■ 10"^° ■ ergs ■ Gauss'^ = 1.16 ■ 10"^^ ■ MeV ■ Gauss'^ . (9.3) 



In one of the traps used in [pSV86|| the magnetic field Btrap is about 20 kGauss and the electric 
field Etrap is about 0.000033 Statvolt/cm=1000 V/meter = 61 Gauss. Then I get 



c • -Pt-rap • J^trap ' ' t- _ ^ ^ i n-21 /ik i n-15 a#„T/ 

9 



7.3 ■ 10"^^ ■ dyne = 4.5 ■ 10"^^ ■ MeV ■ cmT^ . (9.4) 



If the electric and magnetic fields are static and parallel to the spin then the numerical value 
(9.4) leads via (9.2) to the following maximal value of the second order SG force 

2 

^ " , - B - B^-a ^ C4 ■ 1.8 ■ 10-21 ■ dyne = C4 ■ 1.1 ■ 10"^^ ■ MeV ■ cm-^ . (9.5) 
2 • ■ c'^ 

The magnetic field has a gradient Gtrap of about 2.4 Gauss/cm so that 

' ^^''"^ " ^ = 4.5 ■ 10-2° ■ dyne = 2.8 ■ 10"^^ ■ MeV ■ cwT^ . (9.6) 
m ■ c 

For the static, transverse magnetic field when the spin vector a is parallel to the magnetic field 
the numerical value (9.6) leads via (9.1) to the following maximal values of the first order SG 
force: 



e • C2 
2 ■ m ■ c 



■ Vp(a^ ■ Bp) C2 ■ 1.1 ■ 10-20 . dyne = C2 ■ 6.9 • 10"^^ ■ MeV ■ cm'^ . (9.7) 



Thus for the Frenkel case ((5.6) with g ^ 2) the second order SG force is comparable to, but 
still smaller than, the first order force. For the GNR case ((5.8) with g 2) the first order 
force is unchanged but the second order force vanishes. 

For protons the values corresponding to (9. 5), (9. 7) are smaller and the second order SG 
force is much smaller than the first order force. 

My main purpose in presenting these numbers is to compare the different terms of the SG 
force. I do not claim that they give a good representation of the real SG force in a trap because 
my formalism only applies to the semiclassical regime and this may not always be directly 
applicable to traps. 
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10 A transformation of the GNR equations 



10.1 

In this section I deal with transformations which in particular allow the GNR equations to be 
transformed into the reduced Frenkel equations. This would, for example, make it possible to 
solve the GNR equations by symplectic methods. 

10.2 

The general equations (5.5) in the approximation that the second order SG terms are neglected 
read as: 

^ m 4 ■ V 



e • (c2 - Ci - 2) p p p r.PjjP 
+ ^r^2 S^^-U^ -U^ -dxU^oj , (10.1b) 

^ ' Til 



.p ^ e ■ g f p p _ p p 



e-{g-2) f^P 



2 • m 



• ■ K ■ K - SL ■ • • <) . (/^, ^ = 1, 4) 

(10.1c) 



Equations (10.1) are obtained by setting C3 = C4 = C5 = in (5.5). In the case where C2 = g 
equations (10.1) simplify to: 

K = K ' (10-2a) 

'"^'~'~''^-s:.-u:.u:.d:Fi, (10.2b) 



2-m2 

• p ^ e-^ / p P _ P P 



e-i9-2) 
2-m 



■ (-^^L • ■ Ul ■ Ul - Sl ■ F^, . Ul ■ [/;) . (//, = 1, 4) (10.2c) 



For ci = g — 2, equations (10.2) are the GNR equations and for ci — they are the reduced 
Frenkel equations. 

10.3 

I now demonstrate how to relate equations (10.2) corresponding to different values of Ci. I do 

P P 

this by defining a transformation of X and 5* via the rule: 

xl^X^, Sl^^S^^^Sl,, (//,i/ = l,...,4) (10.3a) 
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with 



U,^X,^U^+^—^.S^^-F,^.U^ , (/i=l,...,4) 



;i0.3b) 



where Aci is a real number. One sees that the spin tensor does not change under (10.3), so 
that this transformation only has an effect on the SG force, as seen below. The transformation 
(10.3) is a straightforward generalization of a transformation given in ||Pla66a|| . The constraints 
(3.8) are equivalent to: 



S^,-U, = 0. (/i=l,...,4) 



;i0.4a) 
;i0.4b) 



On introducing the abbreviations 



x„ = (fKx. 



5 ^4:)tJ. , 

d d d d 

Wax^'ax^'ax;)- (/^ = i'-'4) (10.5) 

I define: Q 

E = E,,(f, -2 ■ X4) , B = S,,(f, -I ■ X4) , 

F^{B,-i-E). (10.6) 

Now I come to the main conclusion of this section: if , S'' obey (10.2), then X,U,S 

obey 

X^ = U^, (10.7a) 



m 4 ■ m'' 

e ■ {g - 2-ci + Aci) 
2 ■ 



+ , ■ S,, -U^-Ux- dxF,^ , (10.7b) 



S^-F^x-Ux-U,-S,^-F^x-Ux-UA , (/i,z/ = l,...,4) (10.7c) 



2 ■ m 
2 ■ m 

where all second order SG terms are neglected in (10.7b). In the special case where Ci = g — 2 
equations (10.7) read as: 

X^ = U^, (10.8a) 

Uu = ■ Fnu ■ Uy ■ Sy^ ■ duF^y + Un ' S ' U\ ' d\F^^ 



m 4 ■ \ 



•^^The partial derivatives 9^ always act on functions depending on Xi, X2, X3, X4. 

•^^My notation is chosen so as to indicate that the functional dependence of E on r, —i ■ is the same as the 

dependence of E^^^ on ,i and likewise for B^^^. Note that F depends on Xi, X2, X3, X4 in the same way as 

p p p p p 

F depends on X^ , X2 , X^ , X^ . 
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+?r^ ■ •S'm. -U^-Ux- dxF,^ , (10.8b) 
2 ■ 



S^-F^x-Ux-U,-S,^-F^x-Ux-UA , (/i, = 1, 4) (10.8c) 



2 ■ m 
e • (g - 2) 
2 ■ m 

where all second order SG terms are neglected in (10.8b). Therefore the GNR equations are 
transformed under (10.3) into equations (10.8). In particular, with the choice: Aci = g — 2, one 
has transformed (10.8) into the reduced Frenkel equations. Thus one has transformed the GNR 
equations into the reduced Frenkel equations. Therefore the GNR equations can be solved by 
solving the DK equations and inverting (10.3b) so that one can use symplectic methods ||BHR| . 



For practical applications in accelerator physics it is helpful that (10.1-8) contain the charge e 
only up to first order. 

It follows from the normalization of the spin vector that: 

5^. ■ S^, = . (10.9) 

As in (1.9) this equation is of second order in spin so that it plays no role in this paper. Note 
also that (10.9) is conserved under (10.7c). 



Summary 

I have studied classical spin-orbit systems at first order in spin and, by applying dimensional 
analysis and imposing Poincare covariance, I have found that these spin-orbit systems are 
characterized by five dimensionless parameters ci, ...,05. My axiomatic approach is supported 
by the observation that the most prominent spin-orbit systems, namely those of Frenkel and 
GNR, are special cases of my scheme. 

In this approach, i.e. at first order in spin, the five parameters are purely phenomenological 
and are to be determined by experiment. For example the Frenkel and GNR equations give 
very different SG forces at high energy in proton storage rings. There are also differences for 
high fields in traps. Theory is of little help. For example, as I pointed out in [o] even the Dirac 
equation cannot deliver unambiguous answers. In the three cases mentioned the parameters all 
depend on g. However, protons, for example, have substructure and one should not assume a 
priori that the dependence of the c's on the g's is the same for all particles. 

In this paper I have concentrated on spin 1/2 particles. Nevertheless my results are formu- 
lated classically so that they can be applied to particles of arbitrary spin. 

In addition I have devoted special attention to the DK equations and have found a transfor- 
mation of the particle coordinates which relates these equations to the Frenkel equations. The 
DK equations are therefore (nonmanifestly) Poincare covariant. The new coordinates differ 
from the original coordinates by less than the Compton wave length and the corresponding 
time variables are the same. Thus one concludes that the particles described by both equations 
are effectively indistinguishable. 

As I have just pointed out, different values of the c's, correspond to different systems of 
spin-orbit equations and can lead to very different SG forces. Thus, before proposing techniques 
which rely on SG forces at high energy, one must decide which equations are applicable. Alter- 
natively one can take the view that a measurement of the forces is in itself a way of discovering 
which equations to use. 
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Appendix A 



A.l 

In this Appendix I derive (2.9) from section 1 and subsection 2.1. I introduce the abbreviation: 

-[f--^]-(^^pA4)) , (A.l) 
and conclude from (2. 4), (A.l): 



> / 



2 



■ ■ • (^M A a) • + • • (7m • ^m)' + • (7m • ^m)' ■ ^m) 



7m 7m 7m 7m 

;;;-(^;;;TTF-^"' ""'* ""■"" 



„2 



6 ^ — * 

m (7^ + ly ^ 



p 



m 7p + 1 m 7p + 1 ^ 7 
, 7p - - 1 - 
7p + 1 

El-Vp-s+ — ■ El-Vp- -Vp-Vp 

^^7p + 1 " V7p 2 J 2 7^ + 1 

-[|-^]-(^^pAip) 

-[|-^].(..A4))-^P-|^-^-^t.,^.[,^,4^4] 
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m 7^ + 1 ^ ^ " m 7^ + 1 ^ ^7, 



+ - — • S ' ■Vp-Vp\-Vp 



7p + 1 

El-Vp-s^ — ^ , • • -y^ • • t;^ • t;^ 



e 

+ — 
m 


r 1 9 
7p 


- 2 
2 




2) 7p 




2-m 


7p + 1 


e 


7^ 




m 


7p + l 


^2" 


e 

+ — 

m 


7p 




7p + 1 




e 


/p 




m 


7p + 1 


^2 


e 


7p 




m 


7p + 1 


e 


1 




m 


7p + 1 


, 7p 


• 


1 

7p 


+ 1 ^ 


e 1 

m 7^ 


El-Vp-S 


. e 


e 

+ — 

m 


r 1 ^ 
• — + - 

7p 


-2 
2 


, e • (5 - 


2) 7p 





Ill Z ~r J. 



p + 1^ 



7p + 1 



2 



-l^].[vl.s-vl-Ep-vl.Vp.s^.Ep]. 

/p • (VpABp+Ep) 

. -Fi , e -fp r 1 



V, 



1 ^ 



m 7p + 1 ^ 7 



2-m •7p--gp-^^p--^'-^^p-^P 



2-m 7, + l ^ l^p/^^p^ 



2-m 7p + l ^ ^ 2-m7p 



- — - S ' - - ■Up 



Introducing the abbreviations 

m ^7p 2 ^ ^ 2-m 7p + 1 ^' ^ -f 

+^^^^-^^-.^-4-(-^pA.1-^---.-'t-4--^P 

2-m 7p + l^ ^ ' 2-m7p 



^ 2-m ^ 2-m 7p + 1 ^ ^ 

Ai = - - [^-^ - 7p - - - (Ep + A Bp) + — - - ^p - s 

mi 7o 



^^Note that all A's are first order in spin. 
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2 2 -fp 2 -fp 

A, ^ - Ao = ^ ■ ((^-2) . -^^-s^-v, ■ (vpABp) 



If 

V (fl -/ , 1 -P 

7p + 1 

one sees by (A. 2) that (2.9) is valid if A2 vanishes. Hence the remaining task of Appendix A 
is to show that A2 vanishes. 



+ {g-2).^—^.{sABp)-{g-2).-^.vl-Bp.{VpA^ 
-{g - 2) . . Bl ■ {sA Vp) ■ Vp) , (A.3) 



A.2 

For the case where Vp,Bp are hnearly independent (e.g. nonparallel) , one has the following 3 
linearly independent vectors: 

Vp, Bp,Vp A Bp . 

One concludes from (A.3): 

v^.A,^ . i;t . ( A 4 ) . [ 1^ + .v^.Vp]^0, 

2-m ^ ^ 7p 7p + 1 

5; ■ A2 = , 

[Vp A BpY .A^^^-iVpA Bp^ ■ ((5 - 2) • ■ ■ Vp ■ {Vp A Bp) 



H9-2).^—^.{sABp)-{g-2).-^.vl-Bp.{VpA^ 
7p 7p + 1 



-(g-2).-^.Bl.isAVp).Vp) 

e-{g-2) _ / 7p . . ^ -liff -v ■B'' -B - B'' ■ v ■ ■ B ] 

+ -[vl-s-Bl-Bp-s^-Bp-vl-Bp] 

7p 

^ -v^ - Bp- y - Vp-B^ -s-v^ ■ s- ■ Bp]] = . (A.4) 

— * — * — * 

Because Vp, Bp,Vp A Bp constitute a basis of vectors, one concludes by (A.4) that A2 vanishes 
for the case where Vp, Bp are linearly independent. 

To discuss the case where Vp,Bp are linearly dependent, one first observes by (A.3) that 

— * — * — * 

A2 vanishes, if {Tp = or Bp = 0. It remains to consider the subcase with: Bp — X ■ Vp, where 
A is a constant which balances the dimensions. Then from (A.3) follows 

= . (iZ^ . (.-A - ^ . . . A .^) ^ . (A.5) 

Hence I have shown that A2 vanishes. This completes the proof of (2.9). 
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Appendix B 



B.l 

In this Appendix I derive (2.10) from subsections 1.1-2 and subsections 2.1-2. This task is 
tedious but straightforward. In fact I only have to consider m ■ {'jp ■ Vp)' as determined by the 
Hamiltonian H^^j through the relation: 

d 

m-irip- Vp)' = {m ■ 7^ ■ t/p, i^^J^ + ^(m ■ 7^ ■ Vp) . 

Hence m ■ (7^ ■ Vp)' is a well defined function of r^^, a and the main task is to reexpress 
it as a function of rp, t, Vp, s. 

B.2 

First of all I express v^^ in terms of t, v^^j, a and I abbreviate: 

A3 = {r-„at . ^^j^^ = f^^^ff^ . wj . (B.l) 
To simplify (B.l) I calculate for an arbitrary function /(vf^,^): 

^„(/ • at . (vf^^ A EJ) = ^„(/ . nl ■ {E,, A <?)) = ■ {n,, A EJ ■ + / ■ (E,, A a) , 
^M^M = ^ , = , (B.2) 

from which follows by (B.l): 



m V " ""J„ 2 ■ 

e/^m^^-^ ^-2^1^ 1 - 

cr ' ■ iJ., • vr., r — -cr'-TT,,-!^ — r • is. 



■ — ; --r - TT' ■rJ.-TT.J ■ ; ■ TT ' ■ -D,, 



•^^Here the partial derivative i9/9t acts on functions depending on fj^^,t,Pj^^,a. 

^"The nabla operator Vj,,j always acts on functions depending on , t, p^^ , a and it is the gradient w.r.t. p^^ . 
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«-2 fe 1 

+ — ^ — --(^ "M 



M V " M 

.-■ + .gt.^ -Tft -5 -Tf +- ^ 7ft -5 -a 



if 2^ " " m ■ (if,, + m)2 ^ m ■ (if,, + m) 

+ ^ ■ — ■ ■ ■ V ■ B ■ — ■ a' ■ V ■ ■ B ■ v 

+ - ■ -fF^ ^l, ■ ■ ^) ■ (B.3) 

Also one gets from (1.4),(B.l-2): 

T- = ^m-Jm = {^M, = {^"m, ^^mIm - {^M, ■ W,,},, = V,, - As . (B.4) 



Next I abbreviate 



Ae = m-^l- vl ■ As , (B.5) 



from which follows by (1.3-4): 



dt J,^ 
= ^ ■ (tt,, A 4J + e ■ 4, + A4 , 



J' = — -E^ ■ 7f,, + vr^ ■ A4 = — ■ ■ 7f„ + {Jt ■ A4 , (B.6) 

•J M M ^" 



SO that one gets from (B.1),(B.4): 



J'., ^ r, e _ - . e - 1 



- • Tf + a;, = ^ ■ (Tf A 5, J + ^ ■ ^„ + ^ ■ A4 - ^ ■ . 



A/ 



■ ■ A4 • 7f,, +/\, = — - (Tf,, A 5,,) + _ . E,, + — ■ A 



- ^ ■ ■ ■ ^A. - ■ 4 ■ ^4 ■ tT,, + A 5 . ( B . 7) 



^The partial derivative d/dt in (B.6) acts on functions depending on r^,, t,PK 
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To eliminate ifj^ from the rhs of (B.7) I use (B.4) to calculate: 

■^M "^M -^M 

777^ -* 

= l-7^ + 2-4-A3, (B.8) 



SO that 



,2 



1 — tT^ • tT 



m 



7^ • f 1 - 2 • 7^ • i;t . A3 V 









1 




■Ae)-^ 


1 

M 


1 

M 


2- Ae 

M 




= ^ _ 


As . 






'^M 



K.. - A„ 



(B.9) 



Inserting (B.4), (B.9) into (B.7) yields 



Mm m m jvi" m 

-e-# •K-A3]-[iT^-(-^ + |f)--^]--^-4-A4-^-^ + A5 



jy- \ M MJ ' jy- M M M M ' jy-2 

^^M ^^M ^^M 

M M MM 

-4^-^^-El-v^-< + ^-El-v^-K,-^.vl.K,.v^+A,. (B.IO) 

M MM 

B.3 

Now I begin to evaluate m ■ (7^ • Vp)' and it is clear by (2.3b) that: 

m-{-fp- Vp)' ^ m ■ -fp ■ + m ■ -fl ■ vl ■ ■ Vp . (B.ll) 

The remaining task in this Appendix is to reexpress the rhs of (B.ll) in terms of rp, t, Vp,s in 
order to demonstrate that it is identical with the rhs of (2.10). 
Introducing the abbreviations: 

At=Vp-V^ , 

As = A'7 , (B.12) 
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one gets by (1.3), (1.6), (2.3a), (2.5), (A.l): 



M M 



e -g 1 - 



At 



2-^ 7m 

1 7_ 



^ 7m + 1 



1 

m 

e 

e 



■ (^Ai?^) 
1 



1 



M M ' 

' ' 1 

m2 7^ + 1 
1 



(^Att^) 
1 



m2 7j^^ + 1 



(^ATf^)' 



# 1 

" " (7m + 1)^ 

# 2m 

" " (7m + 1)^ 



(^A7r^) + -^ 



7m + 1 

(^A^;^) + -^ 



m2 7^ + 1 



(^Att^)' 



7a 



M M 



(7m + 1)^ 



m2 7j^^ + 1 



m 

7m 



2 ^L-^M 



1 



(7m + 1)^ 



1 



m2 7„ + 1 



((1^ Aa) A7f^+aA7f' 



M M 

/M 



7m + 1 



« ' • ^M • ^^M. 



• ^M • M + [ • + • • ^^M • ^m ] A ttO 

7m 7m ~r 



1 



■ V 

M M 



{1m + 1)^ 



e • g 1 
2 • m2 7j^^ + 1 

1 1 . 1 



(sA^;^) 



IM IM ' 



-Vm- Vm) 



7a 



M. ^1 .if . Q 

m 7mH-1 " " 



7a 



O , L •S + 

^ 7m + 1 7m 7m + 1 
1 



• ^M • ^m] a [e • (iv A S^) + e • 



m 

e ■ fi- 1 



2--^I-^p- (^^^l)2 -(^A^^) 



2 -m^ 7^ + 1 ^ ^ 



+ 



7p 7p + 1 



s'' -Vp ■ Vp] - 



7 



w? 7p + 1 7, 



7p + 1 



m 7p + 1 



^2 - ^ (7p + l)' 



s A V, 



e-g 



2 • m2 7^ + 1 



7 



e 



7^ 



m? 7p + 1 



s' ■ V, 



If 

K ^ 

•7p 2 J 



P ^1 ^ ^ 1 
s ' • Vp • VpJ 



7p + 1 



^2 

e 



7^ 



7p + 1 



(^pA^;^) 



1 



m2 7^ ■ (7^ + 1) 



2 7^ + 1 ^ 



7p D -*t ^ 



"^■7p-(7p + 1)'^ "^^■(7p + 1) 

e 7p ^+ ^ ^ n 6 1 

— 2 ' 7 — rT\^ • ^^p • f^P • s ^ • • 5^ + — ^ — — 

(7^ + 1)^ ^ ■ ['jp + 1) 



isAE. 
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2-m2 7^ + 1 ^ ^ ^ ^ 2-m2 7^ + 1 " ^ 

+;^-^-^-(?A^„). (B.13) 

Next one concludes from (B.10),(B.12): 



jy^ ' • M / • jy- M jy M M ^ M • t^2 

iA,ABj + 4-.Ae-E^ + ^-A, + ^.El.A,.v, 



-^■^^■El-v^-v^ + ^-El-v^-As-^-vl-A,.v^+A, + As. 

(B.14) 

B.4 

If one inserts (B.14) into (B.ll) then the rhs of (B.ll) depends explicitly on ■y^. In this 
subsection I show how the variable iT^ on the rhs of (B.ll), (B.14) can be replaced by Vp. I 
abbreviate: 

A, = -Kp-jl-vl-Ar. (B.15) 

From this follows by (B.1),(B.12): 



1 



A • (1 - 4 • vj = :i ■ [1 - (v, - A,)t . (vp - A,)] 



M 

= ^,-[l-vl-Vp + 2.vl.A,]^^^.[l-vl.Vp].[l + 2.^l.vl.A,] 

1 r 9 ^+ r n 1 r 2 ■ Ag, 1 2 ■ Ao 

= ^.ll+2•7^.^A,]=^.ll-^] = ^-^, 

p p p p p 

^p 

K^=Kp+Ag, 

< r 1 Ag, ;^ , 1 ^ Ag ^ 



= — -A,. (B.16) 



X„ 'Kp L P 'J ^ ^2 P ^ 

M P p P p P 

Also one has by (2.5), (B.1),(B.3),(B. 12-13), (B.15): 
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p pi 



m-{Kp+m) ^ '7^ 7p + 1 

fi'-2/e 1 em + 2-ir„x^ ^, 
H • — • ■ s' ■ ■ B^, • -, r— ■ ■ v„ ■ ■ B,, ■ v„ 

2 \m Kp + m ^ m {K^+m)^ p p m p 

g — 2 / e 1 ^t^^ e 1 ^+^^+^^ 
2 \m Kp + m -fp + 1 p p m p 



m2 7^ • (7^ + 1) 
A3 + A7 = — ■ — - -Bl-s-Vp 



m2 7^ • (Tp + 1) 2 • 72 



' p 

„2 



(7p + Ij 2 2 ■ 

^ ■ -^L ■ 'f^P ■ 7 ■ (^A iTp) % ■ 7 — ■ BI^ ■ Vp ■ s , 

m2 ^ ^ (7p + 1)^ 7p-(7p + 1) 

A6 = m-7^-^;;-A3. (B.17) 
Combining (B. 16-17) I can replace on the rhs of (B.14) the variable v^^ by Vp-. 



r e ^ e-Ag^ e t*, ^ ,6 e-Ag,-^ 



• • • ^P - ^ • ^P - • A7] • [^^P - Ar] + • Ag • (^7, A BJ 

P P P p 

-^■{AsABj + ^-A,-E^ + ^-A, + ^^-El-As-Vp 
e . ^ ^ e -5+ ^ - 1 



^2 A6-Ei-i;,-i;, + — ■Si-iT.-Aa- — ■i;;-A4-iT, + A5 + A, 



'p p 



-^■{VpABj + ^-E^-^-El-Vp-Vp- ^ • ( A ) 
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P P P 



Ae • # • • + • • • A3 - • ^Tl • A4 • + A5 + Ag , 



(B.18) 

so that (B.ll) reads as: 

m- {jp- Vp)' = m ■ 'jp ■ v' + m ■ 'jl ■ vi ■ v' ■ Vp 



+-^.[A,-A,].[VpAB^+Ej-e-[l + ^l-Vp-v^p]-[{As + Ar)ABj 
J^p 

+e. [El- Vp+^l-Vp. El- Vp.v^^+^l-Vp. El]. [As + Ar] 

+ A4 + Kp ■ [As + Ag] + . 72 . . + Ag] • i;, . (B. 19) 

B.5 

In this subsection I show how the variable on the rhs of (B.19) can be replaced by fp. The 

— * — * 

dependence on comes in only via the field vectors E^^ , and their first derivatives. First 
of all by using (2.1), (2. 5) I abbreviate 

m 7^ + 1 ^ ' 
An = (Alo-V^)^^, 

A12 = (Alo-V^)^^, (B.20) 
from which follows: 

= E^{r^,t) = i^(r-; + Aio,t) = ^^(rv,i) + An = 4 + An , 
= ^m(^"'m, = B^{Tp + Aio, t) = B^{Tp,t) + A12 = 4 + A12 . (B.21) 
Inserting this into (B.19) yields: 

m- {-fp- Vp)' = e- {VpA Bp) + e- Ep + e ■ {Vp A A12) + e • An 

^ V ' ^ V ' 

Lorentz force first order SG terms 

+A4 + Kp ■ [As + As] +Kp-^l-vl- [A5 + As] ■ Vp 

•• V ' 

first order and second order SG terms 

+ ■ [Ae - Ag] ■ [vp ABp + Ep]-e-[l + -fl-Vp- vl] ■ [(A3 + A7) A Bp] 
p 

^ V ' 

second order SG terms 

+e-[El-Vp+^l-Vp-El-Vp-vl+^l-Vp-El]-[As + Ar] . (B.22) 

^ V ' 

second order SG terms 
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Introducing the abbreviation 



Ai3 = e • {Vp A A12) + e ■ An + A4 + i^p ■ [A5 + Ag] + ■ 7' ■ ■ [A5 + Ag] • 
■ [Ae - Ag] • [v^ AB^ + E^]-e-[l + 'yl-v,- v^^] ■ [(A3 + A7) A B^] 

+e.[El-Vp+^l-Vp-El-Vp-v^^+^l-Vp-El].[As + A^] , (B.23) 

one then gets 

m • • Vp)' ^e-{VpABp) + e-Ep+ A13 . (B.24) 

I now have to show that (B.24) is identical with (2.10). Therefore the remaining task of this 
Appendix is to simphfy A13. The first two terms on the rhs of (B.24) constitute the Lorentz 
force whereas the remaining part constitutes the SG force. To disentangle A 13 it is important 
to notice that the SG force occurs in two different forms. The 'first order part' contains the 
field vectors Ep,Bp only linearly; more specifically it is linear in the first derivatives of the 
field vectors. The 'second order part' contains the field vectors Ep, Bp quadratically. Note that 
in the second order part no derivatives of the field vectors occur. In equation (B.22) I have 
indicated which of the two forms of the SG force occurs in a term. 

Accordingly one can split the SG force A13 into a first order part plus a second order part. 
First I abbreviate by using (2.8): 

Am ^ -V,(at . ^ _v^(o; . [± . -I^ . . . v^]) 

Ip 7p + 1 



m 7p 



L.v,(...[,±,i^,,,,_^.,,.,,.,,_l|__>^H.,A^,)]) 



m 7p V'7p 2 ^ ^2 7^ + 

dE 



■ E' A\—-s + 



m • {Kp + m) ^ 7p 7p + 1 
^^'Tp-CTp + I)'"^^'^^'' 

^""-^•"^^•^^^•(^;tif-^^^^"^-^-^- 7p-(71 + i) "'^^ 
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The partial derivative d/dt in (B.25), (B.30) and (B.32) acts on a function depending on fp,t. 
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- 2) 7^ ^ ^ e ■ - 2) 7^ ^ 



2-m? 7^ + 1 ^ " ' 2-m2 7p + 1 " ^ 

Al7 = A4 - Ai4 , 

A18 = A5 - Ai5 , 

Ai9 = As - A16 , (B.25) 
from which follows 

A4 = Ai4 + Ai7 , 
A5 = Ai5 + A18 , 

Ag = A16 + Ai9 . (B.26) 

Note that A14, A15, Aig are linear in the electromagnetic field vectors, whereas A17, Aig, 
Ai9 are quadratic. With (B. 25-26) I can now abbreviate 

A20 = e • A A12) + e • An 

+A14 + ■ [Ai5 + Aie] + . 72 . iJt . [A,5 + Aie] • v, , (B.27a) 

which denotes the first order part of the SG force and 

A21 = Ai7 + ■ [A18 + A19] +K^--fl-vl- [A18 + A19] • 

-t\ p 

+e ■ [El .v^+^l-v^-El-Vp-vl + ^l-Vp- # ] • [A3 + A7] , (B.27b) 
which denotes the second order part of the SG force. Then 

A13 — A20 + A21 . (B.28) 

Thus (B.24) reads as: 

m-i^p-vj = e-{VpABp) + e-Ep+ A20 + A21 . (B.29) 

first order SG terms second order SG terms 



B.6 



In this subsection I simphfy the first order part A20 of the SG force. First of all I calculate by 
using (2.8),(B.20): 

VpAiVpABp) = -{vl-Vp)Bp, 

BE 

• Bp) = (4 . Vp)Bp + Vp A {Vp A Bp) = (4 . V J4 +VpA-g^, 
e-{VpA A12) = e-VpA [{A\o ■ Vp)Bp] = e ■ (AIq ■ Vp){Vp A Bp) 
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= e • Vp[Alo • {v^ A B^)] - e • Aio A [V^ A {v^ A B^)] 

^ ■ V,((5 A v,)^ ■ {v, A 4)) + e • Alo A [{vl ■ W ,)B^ 



m 7p + 1 

m7p + l ^ m7p + l 



m 7^ 



1 / dE 
——^i^-s^.v,-{vl-V,)B^-s^-v,^ 

^ • ^ • [(.^A.jt . v,]4 = ^ • ^ • • (.^A..)) 

(-V,[il • A - (.-A A ^) 



e • An = 



— * — * 

Here I used the fact that the spatial derivatives of Ep and Bp only appears in the SG terms, 
i.e. in leading order spin. Therefore one can always approximate: 

dEp. dE... OB p. dB... 

P,k M,k P,k Al,k / ■ ^ ^ ^\ 

dr dr ' dr dr ' ' / 

Secondly I conclude from (B.25) 

Ai5 + Ai6 = — ■ — - -s' -B -Vp 

Tp ■ (7p + 1) 

+ ^7p-(7p + 1) 2 J ' 



e . 7p + 2-7p ^7. .r^' atT 1 I ^•(^~^) ?t . o/ 

-4 ■ ^P , ii,2 ■ (^"/^^^p) - 4 • r— Tn 

^ ^ {'jp + 1)^ 7p ■ {'jp + l) 

>,2 



X,-[Ai5 + Ai6]+i^p-7p-^p-[Ai5 + Ai6]-^p = 3- -s^-B'^-Vp 



e_ _ 7p + 7p + 1 
m 7p + 1 
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m (7 + 1)2 VP 
+-•[ — + ^\.{E' As) 

Combining (B.25),(B.27),(B.30-31) one gets 



A. 



1 / r)F 



+(5*, • v,)(5', ■ B, ■ si - St ■ v,(5-' ■ B, • 



e ! + - + 1 



m 7p + 1 



IP 

m (7 + 1)^ ^ 
+-•[ ^ + ^_].(^' As) 

+ m ^ (7p + l)2 +2 ^ l^pA^pJ+ ^ - 

--■v^. E' ■ , • (s A - - • -^-B 



' s 



m " ^ (7^ + 1)2 ' m 7^ + 1 ^ ^ 

— * 

• — • V, st • 5J + — - . . s - . ^ . g- s-t . 5^ . ^ 

2-m m7p + lV at p p 

dB„ ^\ e 1 / ^ rf^t - M 9B„ ^ 95„ ^ \ 



„2 



+ — n — --yp 



„3 I o ^,2 



e 7^ + 2-7p -.t 



m (7 + 1)" 



s^-{VpAE')-v^ 



+-•[ ?_ + ^_].(^' As) 

m ^ 7^ + 1 2-7^^ ^ ^ ^ 
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^ (7, + 1)2 +2 ^ l^p/\t;,j+ s V, 

m ^ ^ (7p + 1) m 7p + 1 ^ 

f^■-•V,(st.5J + ^:i^.7,■.-t.^;^.5; + ^.7,■5't.S;.^;, 
2 ■ m '-fp 2 ■ m ^ m 

■ --^As' ■ (v^AEp)] - - ■ + ^ • ■ . ivpAE'^).Vp 
2 ■ m '^p m (7 + ly 

m ^ 7^ + 1 2-7^J ^ ^ ^ 

m ^ (7p + 1) 



This can be further simphfied by calculating 

vl-E'^- (s A Vp) -vl-Vp- (s A = s-A (t;^ A (^7^ A E'^)) 

= Vp-s^ ■ {Vp A E'^) -s^ -Vp- {Vp A , (B.33) 

from which follows 

^2, = ^---Vp{s^-Bp) + '-4^-l,-s^-Vp-B'^ + ^-^p-s^-B'^-Vp 

2 ■ m ^ m ^ 



e 


■ 9 


1 


2 • 


■ m 


7p 


e ' 




1 


2 ■ 


m 


7p 


e ' 


■9 


■7p ■ 


2 • 


m 



. .Vp[s^ -{VpAEp)] 



^pJ ■ 



m ^ 7^ + 1 2-7^ ^ ^ (7^ + 1)2 2 ^ ^ ^ ^ 

^ (7p + l)2 (7p + 1)^ 2 ^"'"''-^ 
^ ■ ^ ■ V, f.-t -Bp- El- {sA Vp)] + 'lp^.(g. [,-t . B' - {sA Vp)^ ■ E' ' 



>j Vp 



2-m-7^ "^V ' ^ ' "^V 2-m 
+ - 2) • + A 5; - • ■ t;,] A . (B.34) 

Inserting (B.34) into (B.29) one has thus obtained: 

m ■ (7^ ■ Vp)' = e-{VpABp) + e-Ep+ ^'^ ■ Vp (s^ • 5^ - # ■ (s A Vp) 

■ ■ ■ " ^ "^-^^ ■ ■ + - 2) ■ [4 + t^p A 5; - ■ 4 ■ Vp] A s 

+A21 . (B.35) 

With (B.35) I have simplified the Lorentz and the first order SG terms and have derived the 
first order terms of (2.10). In the above derivation it was essential that the electromagnetic 

"^^One sees by (B.35) that by neglecting second order SG terms the charge e appears only up to first order. 
Thus one could have derived the first order SG terms in an alternative way by making first order perturbation 
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field obeys (2.8), i.e. is a solution of the vacuum Maxwell equations. The second order SG 
terms are simplified below. 



B.7 

In the remaining subsections of this Appendix I complete the derivation of (2.10) by disentan- 
gling A21, i.e. I have to deal with the second order SG terms. Q 

First of all I simplify the rhs of (B.23) by collecting its terms in a convenient way and to do 
this I calculate by using (2.12),(B.5),(B.12): 



(^K^ ■ [A3 + 7^ ■ tJt . A3 . v^]^ = ^- (^K^-[m- As + m-^l-vl- As-v, 

= — ■ K' ■ [m ■ A3 + m ■ ■ if ■ A3 ■ Vp] + — ■ ■ [m ■ A5 

+m ■ 7^ ■ iTj, ■ A5 ■ + m ■ 7^ ■ ■ (7^ ■ v^)' ■ t;^ + m ■ 7^ ■ tTj, • A3 ■ (7^ ■ Vp)'] 
e-El-Vp-A3 + e--il-El-Vp-vl-A3-Vp+Kp-A^ + Kp--il-vl-A^-Vp 

+e ■ 7^ ■ A^ \vp ABp+Ep]-Vp + e-^l-vl-A3- [v^ A 5^ + E^] , 

Kp ■ [A7 + 7p ■ vi ■ A7 ■ Vp]) = - ■ f ■ [m ■ A7 + m ■ 7^ . {Jt . A7 ■ t;/ 



m 

= — ■ K' ■ [m ■ Aj + m ■ -fl ■ vi ■ A7 ■ Vp] + — ■ Kp ■ [m ■ As 

+m ■ -fl ■ vl ■ As ■ Vp + m ■ -fp ■ A\ ■ {-fp ■ v^)' ■ t;^ + m ■ 7^ ■ {Jj, ■ A7 ■ (7^ ■ Vp)'] 
e ■ El ■ Vp ■ Aj + e ■ ■ El ■ Vp ■ vl ■ Aj ■ Vp + Kp ■ As + Kp ■ ■ vl ■ As ■ Vp 

+e--fl-A\- [Vp ABp+Ep]-Vp + e--fl-vi-A7- [Vp ABp+ Ep] . (B.36) 



Thus (B.23) can be rewritten as: 



Ai3 = e ■ {Vp A A12) + e ■ An + A4 + [Kp ■ [1 + -fl ■ Vp ■ vU ■ [A3 + A 



first order SG terms 



first order and second order SG terms 

- e ■ [A3 + A7] A Bp . (B.37) 

second order SG terms 

To obtain A21, i.e. to identify the second order SG terms on the rhs of (B.37), one first observes 
by (B.25) that A4 contains the second order term A17. Also one has by (B.5),(B.12): 

Kp. [l + ^l- Vp. vl]. [A3 + Ar]J = Kp.[l + jl.Vp.vl]-[A, + As] 
+ [Kp. [1 + ^1. Vp. vl]J .[A3 + A,]. 

On the rhs of this equation only the first part contains first order SG terms and by (B.25) one 
observes that the first order terms are given by Kp ■ [1 + 7^ ■ -1/^ ■ tTt ] ■ [A15 + Aiq]. Introducing 

theory w.r.t. the charge. This approach is chosen in ]DS70| , so that from this point of view the first 6 subsections 
of Appendix B are just a check of [ DS70 |. 

'^''To the knowledge of the author this is the first treatment which takes the second order SG force into 
account. 
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the abbreviations: 

A22= (i^p•[l + 7^^?p•<]•[A3 + A7]y-/^,•[l + 7^^Tp•^7;]•[Al5 + Al6] , 
A23 = -e • [A3 + A7] A , (B.38) 
one thus can simphfy the second order SG terms as follows: 

A21 = Ai7 + A22 + A23 . (B.39) 

B.8 

In this subsection I simplify A22 and I first of all calculate 

^ Tp ■ hp + 1) 

, e 1 g-2 ^ w ^ 



m2 '(7p + 1)' 2-7 

■ I- 41 ^ ^ 

e-^i^ 1 /^^.e 1 ^+^^ 

^ (7^ + ly 7p ■ (7p + 1) ^ 

^P • [1 + 7p • ^P • ^] • [A3 + A7] = — • • • • 



e 7p + 2 ■ 7 ' 
m (7 + 1)^ 



+- • 7 + 77^ • {Ep A s 

m ^ 7^ + 1 2-7^^ ^ ^ ^ 

+ rn ^ (7p + 1)2 2 ^ l^p ^ j + ^.^ 7p « ^^p ^p 

--■v^ - E- , ■ (s A - - • • . s ■ s , (B.40) 

m ^ ^ (7^ + 1)2 ^ m 7^ + 1 ^ ^ ' ^ ' 

from which follows by (B.38): 

^ ^^^^■^^!-7p]-.--.^p-^pyAi.+^.r7p-.--.^py-^ 



p 



m (7p + l)' ' 2 '^^ -^P py P . 

■- • f , • (s A iL) • iJt V • 4 - - • f ■ s- v^]' ■ Bp . (B.41) 

m V(7p + 1)2 ^ ^ m V7p + 1 " ^ ^ 
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To simplify this, I calculate by (2.3a), (2.5), (2.9), (B. 12-13), (B. 18): 



^ ^"^ • (sA Bp + {VpAs)A 4 



A Bp] - vl ■ Ep ■ s - ■ ■ Vp ■ vl 



2 • m • 7p V 



2 • m • 7p , , . 

-El-Vp-s^'-Bp-^l-s^-Vp-El-Vp-v^^-Bp^ 



e-g 



2 • m • 7, 



.Ei-s-v^-Bp + 



•(^-2) 
2-m 



• 7p • I s 



■Ei-Vp- vi ■ Bp\ + ^— • it • iJ, • • 5, , 



-S -Vp-j^ -Up-v j^p I , ^ 

^ ^ ' m • jp ^ 



IP 

{s A Vp)' = (7^ -{^AVpij = (^7^ • {a A vj^ = m • (7^ + 1) • A7 

/ 1 \ / 

-7m-(7m + 1)-('?A' 



= m • (7^ + 1) • A7 + - • 



7m + 1 

m • (7, + 1) • A7 + - ■ • < • 4 • (s A Vp ) 

m 7„ + 1 ^ 



M / 



e 


9 1 


2 • 


m 7p 


e • 


(^-2) 




2 ■ m 


e 

^ — 


• — • s ' 


m 


7p 


e 

+ — 


s 


m 


7p 


it . 

p 


{sAVp)' 


e ■ 


(^-2) 




2 • m 


e 

+— 


• — • s ' 


m 


7p 



'^p 



7, 



'p ^p 



"^-•"^-+ 2.m 



e-ig-2) . ^ - 



"^■Bp-Vp 



^-^f-^.^^.s"^ .Vp-{VpAEp) 

Zi ' Tib 



2-m -fp ^ ^ 



7 



P • • t^P 



^4^-7p--'t-^-ii- 



EpAs'^EpA 



^ ^'^ ■ (sA Bp + {Vp A s) A 4 
2-m-^p \ 



[Ep +VpA Bp] -v^ ■ Ep- s--fl- - Vp- vl 



+1^ 7p -s^-i^P • 

2 • m • 7p V ^ 

^4^-fe-^p-«"'-«"'^-4-4+4-4-(^^pAs-)-ii-(iTpA5-) 

Zi III I p \ 
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(st . = (^t . = . A a) + at . 
= — • (s A v^y -n^ + ^-a^ -[VpAB^+E^-Ei-v^- v^] 

e 1 ^4. ^ ^+ ^ 
■ Vp ■ El ■ Vp 

m ^ 

e-(g-2) 1 ^ - e-(o-2) 1 , ^ ^ 

+--i--f-^l-4--. 

m 7, 2 73 

Sj-Bi-v'p = s,-^^-^^-Bl.Ep-^.^^.v^^.Ep.v^^.Bp) , (j = l,2,3) (B.42) 

from which follows by (B.41),(2.9): 

X 7p • (7p + 1) • (2 • 7p + 1) - it + 7p + 1) • (2 • 7p + 1) ^ r - 

A22 = — -^T^ Vl -Ep-S^ -Bp-TT^ 

e 7p + 7p + 1 R ^ ^ 7p + 7p + 1 ^ 
+;^- 7p-(7p + 1) -^ ''-•"-^m- 7p + 1 

(2-7p + 2)-(7p + ir-2-(7p + l)-(7^ + 2-7p) , ^ . . E • vf 

(7p + 1)^ ^ ^ ^ ^ ^ 

e 7p + 2 • 7p .^.^Nt £^ ^/ e 7p + 2 ■ 7^ ^ -.^ /-a-V 

(7p + 1)^ 2 ■ 7^ ^ m 7p + 1 2 • 7^ 

(7p + 1)^ • (2 • 7p + 2) - 2 ■ (7, + 1) ■ (7^ + 2 • 7,) ^ 

w? ^ (7p + l)2 2^ p M p 
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^ ^ m 7p + 1 ^ ^ 



m 7p + 1 

2-7^ + 1 ^ g - , ^2 7p + 7p + 1 R r- A R . p 1 



_ 7p + 7p + 1 
m 7p + 1 



2 ■ m ■ 7p 



+ 



2-m 

• # • Vp • s t • 5 



^ • {Jp • # • - s ^ • iTp • # • TTp • ij; • 



m ■ 7p 



7p + 2 ■ 7^ 



(7p + 1) 



'p --yp 



(7^ + 1)^ 



m (7p + 1) 
e 



2-m 7p 



■ - 2) ^ ^ 



m ^ J (7p + 1)^ 2-7^ ^ 



^^.(El-Bp.s-s^.Ep.B 
.2 - m - 7p V ^ 



e r 2 ^ 

H \ — 

m 7p + 1 2-7^ 

As) - -fl ■ ■ Vp ■ El ■ Vp ■ {Ep Av^ 



^^.{^l.s^.Vp.Ei-Bp.v 



p p p ^p 



-^l-s^-Vp-El-Vp-Bp-vl-Ep-{E, 



(7p + 1) 



m^ 



(7, + l)' 2 

2-m 7 m, ^ 7 



e ■ (g - 2) 1 ^ 



/p 



m 7p 2 73 



El-s 



2-m \ 2-m 7p 



7 • T?i . ^ P m ^ V. 



1 5 7p-l- 



2-m 7 
2-e^ 7/ 



^2 (7p + 1) 
1 



m^ 

„2 



^ ^ ^ ' m ^7p 2 



7p 



El-s^ 



-^■(:^;^-(^^A^;p)-#-[.^pA5,+4 
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m (7p + l)2 ^ " 



2 • m 7^ 



e • - 2) ^ ^ ^ e ■ - 2) ^ - 



2 • m 
e 1 



2 • m 



m 2 ■ m 

± . _L . fit . e _ £ . _L . . £ 7f . B 

P P P P P P 

m ^ m ^p ' 



p p 
e 1 
m 7p 



m 7p + 1 
m 7p + 1 



e-g 



2 • m • 7^ 



sASp + (i;p As) A£;^ 



e • - 2) 
2 • m • 7p 

1 



v2 . ?t 



■{-fp-s' -Vp- 



[Ep +VpABp]-v^-Ep-s--fl-s^-Vp-v^p-Ep-v, 



m? (7p + 1)' " " 
1 



^Jl • • s t • 5, • + • + - P -Bp- [Vp ABp+E^ 



m2 7p • (7p + 1) 



-g 



■ El ■ s- • B^ ■ 



2 " ^p ^p 



+ 



2-m? (7p + 1) 

w? U 2 (7, + 1)2 7^ + 1^ 



2-7, 



w? (jp + 1)3 
7' +2-7, 



(7p + 1)2 
4 7^ 
4 7"^ 
4 7p 



■ •(sAi;p)t-^p-(^;pASp) 



7p • (7p + 1) 


/p ' - /PI 
(7P + 1P^ ^ 


1 


7p-2 


7p • (7p + 1) 


7p • (7p + I? 


1 





El ■ A Vp) ■ Ep 
■vi-Ep-{sAEp) 



■El-Bp. s 



+— • \-— + 



, g 2.^^ + 2.^^^^,-i 



4 2 



ll ■ (7p + 1) 



-Ep-B^ 



U 72 ^ 7p-(7p + 1) ' (7p + 1)^^ 



+ 



■El-Ep-{VpAs) 



-g 



1 



2.m2 (1 + 7p)^-'''-'^^-''^-'^^-''^ 



2-e2 



7f 



+ T ' "..r i • {sA Vp)^ ■ Bp ■ (Vp A Ep) 



m? ijp + 1)3 

^ .g^_g_ 2-7^+4-7p + l 7^ + 2-7p . 
w?'U 2 (7^ + 1)2 + (7^ + 1)2 ^ ■ 

7p-l 7p + 7p-l , 7p + 2-7 



4 7 



7p ■ (7p + 1) (7p + 1)' 
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El.s-{VpAEp) 



B.9 

In this subsection I simplify A17 and A23 and use: 

3 

n=l 

3 

<^k ■ {t^mj^^m.Jm =e-ak-J2 ^imn ■ , (j, A;, m = 1, 2, 3) 

n=l 

which follows from (1.1),(1.2),(B.21). Thus I get by (1.3),(B.5),(B.25): 

Ai7 = A4 - Ai4 = {^M, • ^m\m - Al4 

e f m g — 2 1 ^i^^t-, ^ 

-to A R ^ 5-2 2- Jj,, +m ^ - a n ^ 

n2 



,2 



4 ■126:- 7:^1 -K^'^^"^^'! 



m2 2 (7, + 1)2 ^ ^ ^ l^pA^.j 
.2 
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2 • 7p 7„ + 1 V ^ 7p 7p + 1 



2 

t 

^2-7, 7^ + 1^ V - " 'i, ' ' ' ' 



Next I calculate by using (B.17),(B.38): 
V Yp- (7p + 1) 



• 1 ■ V _ — I— I I ^ 

2 



Sp-(7p + 1) 2 ^ P P P P 

r 1 9_ 

m2 '47^ + 1)2 2-72 



p 



7p + 7p + 1 



m'^ 7p • (7p + 1) 

-^•[^-^^-^]-^^-^p-l-^^-(-^pA^J 

-^•[(^;il)^-2^]-^^-(-^.A4)-(^^pA4) 

2 



[El-B.-s-Bl-s-E,] 

^ (7p + 1)^ ■ 2^ ^ .p.[4-4-^p-^l-^p-4] 
e2 1 



m? 7p ■ (7^ + 1) 2 • 72 

r + 2 ■ 7p 



P - 1 " 1 1 
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B.IO 

Combining (B.39),(B.43-45) the second order SG terms read as follows: 

A21 = Ai7 + A22 + A23 



p2 . ^ 

.1 I ■ E ■ (sA V )^ ■ E -v 

m2 (7p + 1)3 ^ " ^ " " 



w? (7p + 1) 



+;SJ-'-T+2 . (T, + 1) 

o2 ^2 



I 1 1 

\ 

p • (Tp + 1) (Tp 



1 



2-m2 (1 + 7^) 
2-6^ 7p 



m? (7p + 1)' 

r^' ^ 2-7p+4-7p + l , 7p + 2-7p i . 5 a p ^ 

+^■^4 •^-^■7p-(7p+i)+t^;tif^'^^'"^'^''^^^ 



/p 

o2 



(7p + 1)^ 

^ 1 

2 • m2 ■ (7^ + 1) 



^ -Ei-iVpA B^) ■ (s A t/. 



+ ^.-7:r—^.-^l-E.-^l-B.-s. {BAG) 
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With (B.46) I have an exphcit form of the second order SG terms and my remaining task of 
this Appendix is to reduce them to those in (2.10). I abbreviate the second order SG terms of 
(2.10) by 

A24 = ^ -[-{g- 2f ■ ill ■ ■ {B^ A s) + {g- 2f -BI-B^- {v^ A s) 

-{g - 2f ■ El ■ s- B, + i-g' . vl . + 2 ■ g ■ 1^1 ■ + 4) ■ ■ B, ■ E, 
-{g - 2f • • 4 • {E^ A s) + {g - 2f ■ {v, A 5 . ■ {v, A s) 
+{g - 2f -El-E,- {v, A s) - (/ - 4 • 5) • ■ B, ■ {v, A B,) 
+{g-2)-g-El-B^-s-{g-2)-g-v^^-E^-v^^-B^-s 
+ig-2).g-v^^.E^.s'' ■B^.v^ + {g-2).g-v^^.s-v^^-B^-E^ 
-ig'-A-g)-El-iv^As)-E^ + 2-g-El-iv^As)-iv,AB^)] , (B.47) 

Therefore the remaining task of this Appendix is to show that the rhs of (B.46) equals A24. 
Introducing the abbreviation 

A25 = A21 - A24 = ■ ({g - 2f . s t . 4 . (^T^ A 4) 
+ {g^-2-g)-vl-v,-s^ -B^-E^ 

8 • 7p " 



+[2 • 5 - (^^1^1)2" ] ■ ■ ■ ^ 

+[9 ■ — + 4 • ^ • J , '-l + lf ^ -K-i'^ ^p) ■ Ep 

+b • + 4 • ^ • ^."^ , ^„ , ■,^2 ] ■v^p-Ep-[sAEp) 



1-7^ 


+ 4' 






+ 7p -1 


7^ 




• (7p + 1) 


1-7^ 


+ 4' 


■9- 


ll 


+ 7p-1 


7? 


7p 


■ (7p + 1) 


1-7^ 


+ 2' 


■9- 


' p 


+ 7p - 2 




7p 


■ (7p + 1) 


1-7^ 


+ 4' 


■9- 


7^ 


+ 7p -1 


7^ 


7p 


• (7p + 1) 


1 








■El-v, 


(1 + 7p)^ 



7p • (7p + 1)' 

7p-l 
ll ■ (7p + 1) 



(7p + 1)^ 



-2.g.j:^^^^.s^.v,.El.v,-B,-8-j^^^ 

+[9' -^-9 - ''^1:1 ' ^ ' + ' "(t + 1)-^^ ^ • ^''^ • • ^ 

2 7p-l , 7p+7p-l , 4-7p + 8-7p i - a p ^ 

+ 19 ■ 4 • ^ • — — + . ^ ^ -El-s-iVpA Ep) 

Tp Ip ■ (7p + 1) (7p + 1)^ 

-{g - 2f ■ Bl ■ (5 A i;,) • 4 + 8 • j;^^^ ■ ■ E, ■ vl ■ E, ■ (s A v,) 

+[/ - 4 • ^ + ^ "^^^ + ^^2" \-^\- ^B,). {sA V,) 
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(7p + 1)2 ' p p p - 
+{g - 2f - vl - B^ - {B^ ^s) - {g - 2f ■ Bl - ■ {v^ ^s) 

-{g-2)-g-v^^-E^-s''-B^-v,-{g-2)-g-vl-s-vl-B^-E^~^ , (B.48) 

— * — * — * 

I have to show that A25 vanishes. One can separate A25 into a magnetic part A26 plus an 
electric part A27 plus a mixed part A28- Therefore I abbreviate: 

A25 = A26 + A27 + A28 , (B.49) 

where 



A26 = 



4 • 



((5 - 2)2 . st . . {v^ A 4) -{g- 2)2 • §1 ■ {s A v^) ■ B 

+{g-2Y.if^.B,.iB,As)-{g-2f.Bl.B,.{v,As)^ 

e-'-jg- 2)2 
4 • 

-Bl-B^-{v^As)) , 



t . ■ {v^ A 4) - Bl ■ {sA v^) -B^ + vl-B^- {B^ A 



s 



p p \ p 

A.-^-(-^-1-4-(^^A^^.)^-4-^. 

+ f + 4=- 9- -7T . , .S2 -El-isA Vp) ■ Ep 

Tp Ip ■ (7p + 1) (7p + 1)2 

1-7^ , , 7^+7p-1 4-7^ + 8-72+8.7p-8 ^ ^ . . 

+ f ^ + 4 • y • ^^-^ — — ^ . , -,X2 ■v^^-Ep-[sAEp) 

Tp Ip ■ (7p + 1) 7p • (7p + 1)2 

2 l-7p , , 7p + 7p - 1 4-72 + 8-7^ ^ ^-^ a 

(7p + 1)2 

-8-^:^-^-^-'^-^^p-^;-t?p-(4A^p) 

72 ' 7p-(7p + 1)+ (7p + 1)2 ^ ^ 



+8-^;^-^-.l-^p-.l-4-(«-'AO), 
A28^^-((/-2-^)-^-;-^-,-5t.4.4 



+ [2 • ^ - ■ (^"/^ ^^p)' ■ Ep ■ {Vp A 4) 

+[/-^ + 2.,.i±^].El.4-^+2.,.^:^pl-..-t.4.5; 

7p 7p • (7p + 1) 72 • (7p + 1) 

~2 • ^ • n / ^2 • s ^ • ^P • • ^P • ^P 

(l + 7p)2 
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-(S - 2) . a ■ tit . . st . B, . ff, - (<, - 2) . a ■ iJt . s - 5* . B, . S, ) . (B.50) 



Hence the remaining task of this Appendix is to show that A26, A27, A28 vanish. 



B.ll 

In this subsection I simpUfy A26- I use the same method as in section A. 2. If Vp, Bp are linearly 
independent, I have the following 3 linearly independent vectors: 

— * — * 

Vp,Bp,Vp A Bp . 

One concludes from (B.50): 
< • A26 = , 

• A26 = , 

{Vp A BpY ■ A26 = ■ {Vp A Bp)^ ■ (fit . Bp ■ {Vp A Bp) - B^ ■ {sA Vp) ■ Bp 

+vl ■ Bp ■ (Bp A s) -Bl-Bp- {Vp A s)) = ' ■ [s^ ■ Bp ■ B^ ■ Bp ■ ■ Vp 

-s^-Bp-vl-Bp-vl-Bp + vl-Bp-vl-Bp-s^ -Bp-vl-Bp-Bl-^ 

-Bp-iP -Vp-B^ -s + B^ -Bp-iP ■ s- iP ■Bp]^0, (B.51) 

p P p P p P P P P -PJ ' ^ 

— * — * 

so that A26 vanishes, if Vp^ Bp are linearly independent. 

If Vp,Bp are linearly dependent, then Vp A Bp vanishes so that one gets from (B.50): 

A26 = ■ • Bp . {Bp As) -Bl-Bp. {Vp A 5)) = . (B.52) 

One observes by (B.52) that A26 vanishes if "tTp = or Bp = 0. It remains to consider the 
subcase with: Bp = X ■ Vp, where A is a constant which balances the dimensions. Then the rhs 
of (B.52) vanishes. Therefore A26 vanishes in any case. 



B.12 

— * — * 

In this subsection I simplify A27. If Vp,Ep are linearly independent, one has the following 3 
linearly independent vectors: 

— * — * 

Vp, Ep,Vp A Ep . 
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One concludes from (B.50): 

ll ' 7.-(7. + l) (7. + l)^ ^ ^'"""^^ 

-I 9 9 , -I 



7p 7p • (7p + 1) 



7p-(7p + 1)- ]-;-^.-(^AEjt.,^j 

-(4-7p + 8-7p)-(7p+7p) + (4-7p + 8-7; + 8-7p-8)-(7p + l)] = 0, 
2 l-7p , , 7p + 7p ^ 1 4-7^ + 8 -7^ /-A-^ pt p 

■^^'■'■7p-(7p + 1)~ (7p + 1)^ ]-^.-^.-KA.1 

+8 • • ^; • ip • ^; • 4 • (^"a • 4) - q , 

4-7^ + 8-7^ + 8-7p-8 ^ ^ ^ 

2 l-7p , , 7p + 7p - 1 4 -7^ + 8 -7^ -.^ ^ ^ 

•^^'•^•^-T^- (7. + l)^ ]-4-^.-KA.) 

-8-^:^-^-«-'^-^T,-#-iT,-(4AO 

7^ %p-(7p + 1)^ (7p + 1)^ ^ ^ ^ ^ 

e^^M? '^^^ '^^ '^^^ '^^ ■ ^'"""^^O 

/ 4-7^ + 8-7^ + 8-7,-8 4-7^+8-7, 
4 • m2 ■ V 7p • (7p + 1)' (7p + 1)' 

+8- ,,-.-t..-J.#.4..-t.4.,^t.,+ [4.7^ + 8-7^ + 8-7p-8 

(7p + 1)=^ " " " " " " ^ 7p-(7p + 1)' 

4 ■ 7p + 8 ■ 7p 3 7p ^ , ^ -t F -^ n 

- (7p + i)^ -8-(^;TIF-"--"-^-"--''--^--''--''--V="' 



— * — * 

so that A27 vanishes, if {Tp, Ep are hnearly independent. 



(B.53) 
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If Vp, Ep are linearly dependent, then Vp A Ep vanishes so that one gets from (B.50): 



^27 



A-w? \ 72 7p ■ (7^ + 1) 

4-7^+8.7^ + 8-7,-8 ^ . ^-aP^^F^ ^ " ^ 4 7^ + 7p - 1 

7p ■ (7p + 1)^ 7p 7p ■ (7p + 1) 

4-7p + 8-7pi ^ . , o 7, 



' ^ '^^ -Ei■E^■{v^^s) + S■^^^■v^^■Ep■v^^■Ep■{s^Vp, 



(7p + 1)^ ^ ^ " ^ " ^ (7p + 1)^ 



(B.54) 



One observes by (B.54) that A27 vanishes if iTp = or Ep = 0. It remains to consider the 
subcase with: Ep = X ■ Vp, where A is a constant which balances the dimensions. Then (B.54) 
reads as: 

X e^-^' f^^^^ A 4-7^ + 8-7^ + 8-7p-8 , 4-7^ + 8.7, ^ ^ . 
^-=4:;;^-(^^^-)-l[ 7p • (7p + + (7p + 1)^ ^-'^-'^ 



4-m2 ^ ^ ^ V 7p-(7p + 1)' (7p + 1)' ^ 

Therefore A27 vanishes in any case. 

B.13 

In this subsection I simplify A28. If Vp,Ep are linearly independent, one has the following 3 
linearly independent vectors: 

Vp,Ep,VpAEp . 

One concludes from (B.50): 



vl- A,s = ■ [{g' -2. g)-vl-Vp-P- Bp- El- Vp 

-'^■9-j:^i;^.-El-s-v^p-Bp.v^^.Vp 

+[9'-2-g- ''l^^l-p -s^-Vp.El.Bp.vl.Vp 

+[9'-'-^ + 2-g- ^''''^^\h -El.Bp.s^.Vp 
7p 7p • (7p + 1) 

^^■9-:^^j;^ys^-E.-Bl-Vp 

+[9'-2-g- ^^^'^l '^p - vl -Ep-vl-Bp-p. Vp 

-{g-2)-g.v^^-Ep-s^ ■Bp-vl.Vp-{g-2)-g-v^^-s-v^^-Bp.El.Vp^^0, 
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+ [2 • ^ - ^"(?^'^ + i)7n • ■ 4 • ^ B,)^ . E, 

-2.g.j^^^^.s^.v,.El.v,.El.B, 
+[/-2-^-i^]-l-4-l-4-^l--' 

-(^ - 2) • ^ • • 4 • s-t . 5, . . _ _ 2) . ^ . ^;t . 5 • ^Jt, • 5, • # . ^ 

(^2 _ 2 . ^) . . t;^ . 5 1 .B^-El-E^ + [g''-2-g]-s'^ -v.-El-B^-vl-E 



4 • m 



^[g^-2.g]^-^.El.B,.El.s + [g'-2.g].El.{v,^B,).{s^v,)^■E, 

^ p 

+[^' - 2 ■ (?] ■ ■ ■ tj; ■ 4 ■ ■ s - - 2) • ^ • ^t; • 4 • • 5^ • ^; • 4 

-(^-2)-5-^^^-s-i)J,-S,-# . (B.55) 



To show that E'J, • A28 vanishes, I introduce the abbreviation 
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A29 = ^b[EI ■ A28] = • ((/ - 2 • ^) • ^Tj, • • # • 4 • s 

+[^2 -2-g]-s^ ■v,-vl-E^-E^ + [g^-2-g]- • # ■ s • 

^g^-2.g].{s^v,)^■E,.{E,^v,)-r[g^-2.g].vl■E,.El.s■v, 
-{g-2)-g-v^^-E^-v^^-E^-s-{g-2)-g-v^^-s-El-E^.v,^ , (B.56) 

so that one gets 



T|, . A20 = • ((.9' - 2 • g) -vl-v.-El-E^-vl 



^^The nabla operator Vb always acts on functions depending on t, Bp, Ep,Vp,s and denotes the gradient 



w.r.t. to Bp. 
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^\g^-2-g\.^^.E^-E\-8-^^-v,-{g-2)-g-^^-E^-^^-E^-^^-s 
-{g-2)-g-v\-s-E\-E^- v\ ■ = , 

+ [g'-2-g]-s^-v^-vl-E^-El-E^ + [g'-2-g]- • # • s • # • ^ 

-{g-2)-g-vl-s-El-E^-El-v\=0, 



{v,AE,)^-A,,^-^.{v,AE,)^.((g'-2.g).v^^.v,.El.E,.s 

+[g^ -2-g]- {sA v,)^ ■ E^ ■ A v,) - {g - 2) ■ g ■ ■ E^ ■ vl ■ E^ ■ ^ . (B.57) 

Hence A29 vanishes, if Vp,Ep are linearly independent. Combining this with (B.56) and using 

— * — * 

the fact that A28 is hnear in Bp, one concludes: 

• A28 = , (B.58) 

— * 

if Vp, Ep are linearly independent. 
Next I calculate 

{VpAEpY-A^s = ^ • {v,/\Epy .(\2.g 
4 ■ 7' + 8 ■ 7p 



ilp + 1)^ 



p ___iLl . (g^ ^Jt . Ep . {vp A Bp) 



+[g^-'-^ + 2.g. ^'^y^-' ].El.Bp.s + 2.g. - ^ . ■ .^t . ^ . ^ 

7p 7p • (7p + 1) 7p • ilp + 1) 



' p 

1 



+ [/ - 4 • g + ^ ] • • (^^p ^ Bp) . {sA Vp) 



7p + 2 • 7. 

(7p + 1)^ 

4 • V 7p 7p • (7p + 1) 



,[^2 o , 2-7^+4-7p + l 4-7^ + 8-7. ^ ^ A ^ )t . 
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+ [-,2 + 4.9-2.9. ^^-i^ - " -^^^ + '^-7- 1 . gt . . gt . flt . (g^ ^ 4) 

-iT!. . . i>; . ij) . (B.59) 

To show that {Vp A EpY ■ A28 vanishes, I introduce the abbreviation 
A30 - V.[(.-. A ijt . A,,] = . ([,^ . iZ^ + 2 . ^ 

2 o 2-7^+4-7, + l 4-7^+8-7, ^ - p 



(B.60) 



so that one gets 



^^•^-^ITT^^At^ •^ + ^-^-^^^- + ^-^- 7p-(7p + 1) 

2 o 2 • 7^ + 4 • 7^ + 1 4-7^ + 8-7^ P A p - -t - ^ n 

-2-^ CrTTip — + {^l + iy ]-<-E,-{vpAEpy-s-vi-Vpj=o, 

^^4o = v^■([^?^^^ + 2■^■^r^^r,+2.,• ^^+^-~f 
4 ■ V 7; ^ 7p • (7p + 1) 

,[^2 o „ 2-7^+4.7p + 1 4-7^ + 8.7p ^ g A ^ )t . . ^ 

+[^ ^ ^ (7p + 1)^ (7p + 1)' ^ " " ^ - " 
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2 I /I . ^, _L 1 /I . ^,2 



2 • 7^+4-7^ + 1 4 • 7^ + 8 • 7 



2 v[2-^- ,2 L 



4-m2 " 7^(7p + l) 

+ [-2-.-^j^ + 4.,..l... 

2 -7^ + 4-7^ + 1 ^ ^ ^ 



+ [-4-. + 2.,.^^J-^ + 2.,.^^^^^ 

= . (B.61) 

Hence A30 vanishes, if Vp,Ep are linearly independent. Combining this with (B.60) and using 
the fact that A28 is linear in Bp, one concludes: 

(VpAEpy-A28^0, (B.62) 

— * — * 

if Vp,Ep are linearly independent. Collecting (B. 55), (B. 58), (B.62) one concludes that A28 
vanishes, if Vp, Bp are linearly independent. 

To discuss the case where Vp,Ep are hnearly dependent, one first observes by (B.50) that 

— * — # — * 

A28 vanishes \i Vp = Q or Bp = 0. It remains to consider the subcase with: Bp = X- Vp, where 
A is a constant which balances the dimensions. Then one concludes from (B.50): 

^2s - ■ {{9' - 2- 9) - vi - Vp ■ ■ Bp - Vp 

~2 • 5 ■ ^ -vl- s-vl-Bp-Vp + [g -2-g- ^ ] ■ s'' ■ Vp ■ iPp ■ Bp ■ Vp 

,[.2 1-7^ ^ 7^+7p-2 ^ ^ 5 I 2 q ~ ^ 5t ^ .B 

ll ^ ' 7p-(7p + 1)^ ^ ^ ^ ' + " " 
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Hence A28 vanishes, if Vp,Ep are hnearly dependent. From this it follows that A28 vanishes in 
any case. 



B.14 



In subsections B. 11-13 I have shown that A26, A27, A28 vanish so that by (B.49) it follows that: 

A25 = . (B.64) 



Thus one has by (B.48): 



24 ■ 



(B.65) 



Inserting (B.47),(B.65) into (B.35) results in 



m- {jp ■ Vp)' = e - {Vp A Bp) + e - Ep + 



e-g 



■Vp[s^ -Bp-El-isAVp) 



2 • m • 7 



p 




2lL.[2-s^ -B'^-Vp-g-isAVp)^ ■E'^-Vp + {g-2)-{E'^As) 
m 

+ {g-2)-vl-s-B'p-{g-2)-vl-E'^- {Vp A s)] 

+ j^-[-(9-'^f-vi-B,-{BpAs} + (g-2f.Bl.Bp.(vpAs} 

^ * fit 

-{g - 2)2 . El ■ s- Bp + {-g^ . . + 2 ■ g ■ v^^ ■ Vp + A) ■ s"^ ■ Bp ■ Ep 

-{g - 2f -vl-Ep- {Ep As) + {g- 2f ■ (vp A BpY ■ Ep ■ {vp A s) 

+ {g - 2f -El-Ep- {Vp A s) -{g'^-A-g)- ■ Bp ■ {Vp A Bp) 

+ ig-2)-g-El-Bp-s-{g-2)-g-vl-Ep-vl-Bp-s 

+ {g-2)-g- .7t ■Ep-s^-Bp-Vp + {g-2)-g-vl-s-v^^-Bp-Ep 

-{g^ - 4 - g) ■ El ■ {Vp As) ■ Ep + 2 ■ g ■ El ■ {Vp A s) ■ {Vp A Bp)]. (B.66) 



This completes the proof of (2.10). 

Appendix C 

In this Appendix I derive the equation of motion (8.15) for m-'-y^j -Vj^j by using those equations in 
sections 1,2,5 and Appendix B which are valid for arbitrary values of Ci, ...,05. 1 only consider 
the case of static (i.e. time independent) magnetic fields and vanishing electric fields. I do 
this for the general case, i.e. for arbitrary values of Ci, ...,05. First of all one concludes from 



(B.12),(B. 15-16): 



m • 7^ • t;. 



m ■ 'yp ■ Vp - Kp ■ {A7 + 'fl ■ Vp ■ V 



(C.l) 
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By (B.13) one has 

At + 7p • • < • At = - T^— TT 'BI'V^-s 

^ ^ 2 • 7p • (7p + 1) ^ 

e-{g-2) 7^ ^| ^ - e 7^ - ^ 

2 ■ 7p + 1 jp + 1 

2 -m^ 7p + 1 
from which follows by (2.9),(B.42): 

(A. + ,1 . . .1 . A.)' = . . 01 . . .-)' 

^' .(.-t.„-,.B,)' + :^,.^.(.-t.B,..-,)' 



2 ■ 7p + 1 
e- g 1 



2-m2 7^.(7^ + 1) 
e • (<7 - 2) 7, 



2 • 7p + 1 



.{if^.B'^.s + Bl-Vp.s') 

(^{s^-Vpy-Bp + s^-Vp-B'^ 



+ —2-^^-{s^-B'^-Vp+Bl.s'-Vp + s^-Bp.^^) 
7p + 1 



2 ■ 7+1 



2-m2 7p-(7p + 1) ^ " 2.m2 7, + ! 



l2 _1_ 1 P P m2 _1_ 1 ^ p ' " P ^pJ 

e-9 7p 



2-m^ 7^ + 1 ^ ^ m2 7^ + 1 



2-m2 7^ + 1 



(^t^p ■ ■ s'^ • t^p ■ Wp + Bl ■ Vp ■ (s^ ■ Vp)' ■ Vp^ 



IP 

w-:^-"^-^--^ •"-•^"-^-^-^-4:^^ 7^(7p + l) •'-•^-'•^'''^-'^ 



2-m2 7^.(7^ + l) " " 4-m3 7^ + 1 

+^^^^-^^-.-'t-i;p-5' +^.^l^.s^.B' -Vp 
2-m2 7^ + 1 ^ ^ m2 7^ + 1 P " 

t?t . 5, ■ . . ;\ 5^ ) _ _y . -vl-Bp-isABp) 



4 • m3 7^ + 1 ^ ^ ^ ' ^ 4 • m3 72 . (7^ + 1) ^' 

One thus has: 

(ir, • ( At + 7^ t/p ■ < • At)) ' = i^p ■ (a^ + 7^ ^P • < • At) ' 
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2-m 7^ + 1 ^ ^ 4-m2 7p + 1 ^ 

2-m 7p + l ^^m7p + l 

m2 7p + 1 p V p 2 • m 7p + 1 ^ ^ 

4 



-P4 ■ ^ZT -4 • -^"p • (^P a4) - • r— TT -^I -4 • (^^ABp) . 

4 • 7p + 1 ^ 4 • 7p • (7p + 1) ^ 

(C.2) 

Combining this with (5.12b) yields by (C.l): 

^ • (Tm • vj' = m ■ (7p ■ v^y - ( ■ (A7 + 7^ . Up . . A7 



= e-(UpA4) + -^-^-Vp(.-'t.4) 

+P^-[(ci + 2)-5-'t-£;-t?p + (c2-Ci-2)-ut .5-5;] 

Z * TDj 

+ 17^ • [C5 • • 4 • (4 A ^ - C5 • S t • 4 • (^P A ^ 

+ (2 • C4 - C3) • S ^ • 4 • (Up A 5p)] + • • Uj, • • S 



4 ■ 7p + 1 ' 



/p 

e • (^ - 2) 7p - R' ^ 7p R' - 

• ■ S ' • Up • 23 • ■ ■ B - Up 

2-m 7p + 1 ^^m7p + l 

2 2 

■ ■ Bp ■ [Vp A Bp) + ^ ■ u' ■ 5 • s ' • Up • Up 



m2 7p + 1 ^ ^ ^ 2 • m 7p + 1 

-'-P^-^-vl-B^-is^v.)^-Bp.Vp 



+f^-^^-^l-4-^"'^-^p-(^pA4) + f^- 7^— TT-^1-4-(«"'a4 

4 • 7p + 1 ^ 4 • 7p ■ (7p + 1) ^ 



e • (Up A 5p) + ^ ■ Vp(st . Bp) 



e ■ C2 
2 • m • 7p 

' s^.B'^.iJ^.ill^ + c^.^p) 

■y„ + 1. 



'2 


• m 




e 


+ 2 


• m 


e 




+2 


• m 


+- 


e2 



^^t-«"'-4-7p-(c2-Ci-^-^.^) + f^-^.ut -s;..- 

7p + l 7p + l 2-m7p + l^^ 



p 

^,2 



• u'f • 5' • s ^ • 
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e2 



■c^-Bi-B^- {v^ A s) 



4 • 

+4^ ■ ■ ^' ■ ^ ■ (2 ■ - ■=» - 

4 • 7p + 1 ' 

Prom (2.3), (C.3) follows: 

"^•(7M-^^M)' = e-(t?p A5p) 

7p 7p + 1 

[^^^ ■ ^ ~ r^XT • • "^p • "^p]^ -B'p-Vp- (;^-^ + ci • 7p) 
7p + 1 7p + 1 

- a-B' --^^ -{02-01 ^ g ■ 

p 7^ + 1 ^ 7p + r 



2 


• m • 




e 


+ 2 


• m 




e 


+2 


• m 


e 


■9 


+ 2 


• m 


e 


■9 


+ 2 


■ m 










+ 4 


• w? 










4 


• m? 










+ 4 





■ • 51 ■ h/p ■ I? ^^^^—r ■ -Vp ■ Vp\ 

7p + 1 



p p 

/p T- -L 

— — ■ v' B ■a^-v„-v„ 

,, I 1 p p p p 



•<-4-(c5-y^- -)•4A[7p•^-^tT•^^•^- 
7p ■ (7p + 1) 7p + 1 

• C5 • • Sp • 7p • {Vp A a) 

•[7p-a--ir-at-^p-t/p]^-5p-(^-pA5j-(2-C4-C3-;i^) 

7p ~r 7p ~r 

\2 „,2 



4 • w? 7p + 1 



IP 

^ ^ %■ Bp - [a AVpY -Bp -v^ 

p 



4 ■ 7p + 1 



,2 „2 



4: ■ m? 7p + 



e -r 7p 
^ ■ 7p + ■ 

e-(VpA Bp) + ^ • V^(a^ • - ^ • • . ^ 5/ 

V p 2-m 2-m 7^ + 1 ^ 



^2.m 7m + 1 7m + 1 

-a-B' ■'y^ -(02-01 9 ■ 

2-m - M /M I 2 1 ^ 

+ ^"^ if! -B' -h -a -a^ -v -v ] 

7m + 1 ^ ^ 7m + 1 
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2-m 7^ + 1 



'4 




4 


• 










+4 


• m? 



■ V 



1 



C5 • • • 7m • (^M A <t) 



[^m-^--^-^^-^^m-^^m]^-^m-(^^mA^m)-(2-C4-C3-;^-^) 



4 •7a 



7m + 



7m + 



4 • 7m + 1 

e'-y-(y-2) 7m 



4 • 



7m + 1 



Ml M M 



.2 ^2 



e -fi- 7 



^•^t-^M-^^-^^M-(^^MA^^) 



4 • m2 7^ + 1 
e.{v,AB,) + ^.V^{a^ -B^) 

Zi ' Tit 



2 


• m 




e 


' 2 


• m 




e 


' 2 


■ m 


e 


■9 


' 2 


■ m 










' 4 


• m? 













B' 

M 


' 1m 


I ' 
7m 


a'^-B' 

m 






7' 

' M 


IM 


+ 1 




■ 

M 


B' ■ 

M 




1m 


M 


■B' 

M 


■ a 


7m + 1 



{c,-g)-c^-—^] 

IM ~ 



+ Ci-7m) 

7' 



7 



7m + 1 7m + 1 



[-^•(^-2)-ci.7m] 



{Bm a a) • 7^ • (C5 - 
■^l-B^-a^ -v^-iv^hB^) 



1m ■ (7m + 1) 



) 



7 



7m + 1 



(C3 - 2 • C4 + Cg) + g"^ 



7 

' A, 



(7m + 1)^ 



+4 



7 



(7m + 1)' 



• C5 ■ ■ ^M ■ 7m • (^V A a) 

4-7, 



7m 



4 ■ rri 



4 ■ m? 7^^ + 1 
e'-g-{g-2) 7', 



Y-^l-5M-(^A^;^)t.S^.i;, 



4 • rm? 7^, + 

This can be simplified by calculating via (1.6-7), (B.4),(B.7),(B. 12-13): 

1 7m 



(C.4) 



^ 7m + 1 



^ m 7^ + 1 



. ^ 1 7 



J. — * J. — * o ~* + 

'"' IM ' ^ ^^M 



Aa) Av^ + a Av 
• -a- B^) 

M Ml 



M 



e 
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_(o _ 2) . ^'^ ■a'' -v^-v^ ■B^-v^ + — -B^ ■ a ■ - — ■ v'' - a- B^] 

\i> ' ^ M M M ' M M " M Ml 

I M ^ IM I M ' 

e 1 / ^ 

" ^ ' ^^~+T ■ ■ ■ • <?+ (5 - 2) • 7m • • ^ • 

from which follows by (B.21) 

m2 ■ 7^ + 1 A 2 ■ " 



e•(^^MA^^)-e•(^;^AAl2) + ^•;^-^•(-|•^T^^^•a 

^^-^-a-^i-^^-^^ + at.^^.^ (C.6) 



9 mMM 

This can be further simplified by calculating via (B.20),(B.30): 

-e • {v^ A A12) = -e ■ [v^ A • V^BJ = -e ■ A^ • V^iv^ A BJ 
= -e • (aIo • (i;^ a S^)) + e • A^o A (v^ A {v^ A BJ^ 



m 7m + 1 "V^ ^ " "V m 7m + 1 

V (^^■•{T -v^ -B -a^-B -v^ - v ) 



I M 



m 7m + 1 

+ - • . [at -B' -v^-v^ -B' -a] 



^ 7m + 1 

e 7a 



m 7m + 1 
so that (C.6) reads as: 



\-v^ ■ ■ V„(at • BJ + a'' -v^-B' -B' ■ - ■ B' -a], 

I M M ' M\ M> ' M M M M M J ' 



e-(iT,A5j = e-(i;^AS^) + ^--^-[at.^^.^^+^t.5^.^^ 

' 7m "T 

-v'^ .B' -a-v^ -v ■ V fa^ ■ 5 )1 + ^- ( -- ■ ■ B -a 

Inserting this into (C.4) yields: 

m-{^^-vJ = e-{v^ABj + ^.-^.[a^.v^.B'^^ 

-4-^L-^-4--^-Vm(^^-^J] + ^-;^-(-|-4 -S^-a 
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' 2 


• m 






2 


■ m 




e 


+ 2 


■ m 


e 


■9 


+ „ 
2 


■ m 








e2 


+ 4 


• w? 










+ 4 


■ w? 


+4 


(7m 






+- 





7m + 1 7m + 1 



I M 
„2 



2-7^ 

^2 

M ^M ""m ^, _L Y ■ [7 I ^ ■ (5 ~ 2) — Cl • 7jj^j 



/M ~r 



7m + 1 7m + 



-B' -a 

M M 



•^L-^M-(^M Aa)-7M-(C5-/ 



1 



7m • (7m + 1) ' 

2 2 

^1L-^m-^^-^^m-(^mA^^)-[-^-(c3-2-C4 + C5)+/- ^r1^2 

7m v7m 

• C5 • • ^M • 7m • (^M A 0^) 



'm 1 _ ^ 



• • ^M • (^^M A ^m) • 7m • (2 • C4 - C3 - ;^^) 

7m 



I • (g - 2)' 7m r^AiT It r r 



" •^t-^M-(^A^-'^)t.i?^,.^- 



2 ^/ _l_ 1 ' 1 "^M/ M M 



4 • 7m + 1 

= e-(t;^A5^) + -^-(-^^-t;i-t;^ + c2)-v^(at.E, 

^ ■ III )J^-\-l. 



[(c2-^)--^-Ci-7m]-^^-^^m-^L 



M 

2 



2 ■ m V 7m + 

+ (ci-7^,+2-7m)-^^-^m-^m 

+-%-[-^-(^-2)-ci-7m]-^^--^m-4-^L--^m 

/M ~r /M ~r 

+(0 - 2) • • • -a 

/M ^ 



e 

4 



M 

2 



4 • rri 



2 



([4-4-7^+7^-(2-C4-C3)]-^J^-a-(^^AS, 



+r^-[-2-C4 + C3 + C5 + 4 + (/-2-^)--^]-4-a-4-S^-(i;^A5^) 

+[(^-2)-^-:r^-7M-c5]-4-s^-(^ASj-c5-7M-^l-^M-(^MAa) 

7m + 1 
1' 



^{g-2f.-^.{aAvJ-B^.B, 



7m + 
7^ 



-{g^-2-g).-^-{aAvJ-B^-vl-B^-v^]. (C.7) 
This completes the derivation of (8.15). 
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Appendix D 
D.l 

In this Appendix I show that (5.12) is equivalent to (5.5). First of all one observes by (3.1-2) 
that (5.12a) is equivalent to (5.5a). 

D.2 

In this subsection I show that (5.12b) is equivalent to (5.5b) and to do that I only have to show 
that the spatial part of (5.5b) is equivalent to (5.12b). To come to that I have to calculate 
the spatial parts of several 4- vectors. Given arbitrary antisymmetric tensors N, N of rank 2 
and a 4-vector a with the following notation: 

— * — * ^ — * — # , 

N ^ {b,d) , N<-^{b,d), a^, = {a\ai)^, (/i = l,...,4) 

one gets: 

N^^-a^= {\aAh + ai-dl\\-a^ -(tj , N^^ ■ N^^ = -2 ■ ■ b ~ 2 ■ ■ d . (/i = l,...,4) 

(D.l) 

One thus gets Q 

P P P 

^tiu ■ Kp ■ Up 

= 7p ■ {[{Vp ABp+ Ep) As~vl-Ep- {v^ A -i ■ {v^ A + E^^ ■ (v^ A s)^ , 
Ux ■ dl = 7p ■ {vl ■ V, + d/dt) , (/X = 1, 4) 

and 

■ ^ 7p ■ [Vp AB^+ E^] , 
SL ■ - -2 ■ V, (st -B.-El- {sA v,)^ , 

K ■ • i^l - -2 ■ 7p ■ [s^ ■B,-El.{sAv,)]-v,, 

K ■ ^1 ■Ul-dlFl^-2.^l-[s^.B'^-{sAv,)^-E';\-v,, 

Sl-K-Ul-dlF:^^^l.[E'^ +v,AB'^-v^^-E'^-v,]As 

"^^The temporal part of equation (5.5b) follows from the spatial part of equation (5.5b) by using the constraints 
(3.8). 

*^The partial derivative d/dt in this equation acts on functions depending on fp^t. 
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■ fZ ■ ■ t/; ^ 7p ■ ( (t^P a 4 + 4 ) a 5, + ■ ■ ) a 

-7^ • {v^ AB^ + ■ E^ ■ {v^ A s) = 7^ ■ [vl ■ ■ {B^ Asj-B^-B^- {v^ A s) 
+El.s-B^-Bl-s-E^+^^^-E^-{E^As)~{v,AB^+E^)^ -E^-iv^As)], 

F;^-SI-F'^^.U; ^^,-[{v,AB, + E,)As-vl.E,.{v,A^]AB, 

-7^ • {Vp ABp + E^y ■{vpAs)-Ep^-fp- [-s^ ■ B^ ■ {v^ ABp) + El-Bp-s 
-(1 + if^-v^).Bl.s-E^-if^-E^-if^-B^- s + if^-E^-Bl-s- v, 
+v^^ -s-vl-B^-E^-El- {v^ A s) ■ 4] , 

F^^-U;; ■Sl-F^,^-2.^,.[s^.B,-El.{sAv,)]-[v,AB, + E,] , (D.2) 

where the expressions on the rhs of the arrows denote the corresponding spatial parts. With 
(D.2) the spatial part of the rhs of (5.5b) reads as: 

^ • 7p ■ (^P A 4) + ^ ■ 7p ■ + ■ V, (i-t . 4 - . (s A i/j) 

2 

+^-[{c, + 2).s''-B'^-v,-C2-{sAv,y-E'^-v, + {c2-c,-2).{E'^A^ 

+ (C2 -c,-2)-vl- s- - (C2 - ci - 2) ■ ■ 4 ■ {v, A s}] 

+ ■ ■ [''^ ■vl-B^-{B^As)-c^-Bl-B^-{v^As) + c^-El-s-B^ 

+ (-C3 • . _ C5 + 2 • C4 - C3) • • • + C5 • i/J, • • (^^ A S) 

-C5 • {Vp A B^y ■E^-{v^As)-c^-El-E^- {v^ A s) 

+ (2 • C4 - C3) -s^ -B^- {v^ A B^) + cs-El-B^-s-C3-vl-E^-vl-B^-s 

+C3 -v^^-E^- ■ B^ ■ + C3 ■ vl ■ s- 4 • ■ 

+ (2 • C4 - C3) ■El-{v^As)-E^ + 2-Ci-El- {v^ A s) ■ {v^ A B^)] , (D.3) 
and the spatial part of the Ihs of (5.5b) reads as: 

7p-(7p-^p)'- (D.4) 

By (D. 3-4), (5. 5b), (5. 12b) one sees that multiplying the spatial part of (5.5b) by m/^p results 
in (5.12b) so that the spatial part of (5.5b) is equivalent to (5.12b). Thus I have shown that 
(5.5b) is equivalent to (5.12b). 



D.3 

In this subsection 1 show that the spatial part of (5.5c) is equivalent to (5.12c). To come to 
that I have to calculate the spatial parts of several antisymmetric tensors of rank 2. Using 
again the notation 

a^ = (a^a4)^, (// = 1,...,4) 
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I calculate first of all: 

1 3 3 

I P P P P P P 

2 ■ Zl ^jki ■ [Skuj ■ Kx -ax-tti- ■ F^^ • oa • Ofe] = 21 ^jki ■ S^.^ ■ F^^ -ax-ai 

k,i=l k,i=l 

3 3 p p p p 

= ^jki ■ [ X! ^km ■ ^mX ' ^X ' Cli + Sj.^ ■ F^^^ ■ Qx ■ Oj] 

k,i=l m=l 

3 3 ^ P P 

— X] ^jki ■ [ X] ^kmn ' • -^mA • ' + Sf.^^ • F^^^ ■ ttx ■ Oj] 
k,i=l m,n=l 

3.3 3 ^ ^ p p \ 

— ^jki • ( 5Z ^^"1" ■ '^n ■ [X] -^mr ' • tti + F^^ • 04 • Oj] + 5';i.4 • ^4;^ ■ ax • aA 

k,i=l m.n=l r=l 

3/3 3 ^ 3^^N 

— 5^ ^jik ' ( 5^ ^nmfe ' •^n ' [ ^mrg ' -^p,g ' Or ' + -^rn4 ' ^4 ' ^i] ~ • F^^ ■ ■ j 

fe,i=l m,n=l 9)''=1 m=l 



3.3 3 

^ ^ ■ ( ^ ^ ^nmk ' ' [ ^ ' ^mrq ' ^ p^q ' (^r ' (^i ^ ' -^p^rn ' (^4 ' (^i 

k,i=l ^m,n=l g,r=l 



3 X 

E Ik ■ ^P,. 

m=l ^ 



m= 

3 3 

^ ! ■ 6ifn Sjm ' ^in] ' ' [ ^ ] ^mrq ' ^ p,q ' Ojr ' ^ ' -^p^rn ' (^4 ' (^i 

i,m,n=l 9>''=1 



3 3 

^ ] ^irg ■ ^ p,q ' Oj^ • (li ■ Sj ^ ^ ^jrq ' ^ p,q ' ' ' 

i,q,r=l n,q,r=l 
3 

-i • 04 • E[-^p,i ■ ' ~ ^p,j ' ' + ■ Ep ■ {qA a)j 

i=l 

— a' ■ s- {Bp A a)j — i • a4 ■ a* ■ Ep • Sj + i ■ a4 • dt' ■ s ■ Ep . + a' ■ Ep ■ (q A a)j , 

(j = l,2,3) (D.5) 

from which follows 

p p p p 

^tiuj ■ -^wA • OA • Oi/ — • F^^X ' (^X ' ^ 

0* ■ s- {Bp A a) — i ■ a4 ■ a' ■ Ep ■ s + i ■ a4 ■ a' ■ s ■ Ep + a' • Ep ■ {q A a) , (D-6) 

where the expression on the rhs of the arrow denotes the corresponding spatial part. For the 
special case: a — U one gets from (D.6): 

P P P P P P P P 

^ixoj • Kx -Ux -U^ - S^^ ■ F^^ -U^ -U^ ^ 

-7p -s^-Vp- [Ep +VpA Bp] + vl ■ Ep ■ s + ■ ■ Vp ■ vl ■ Ep ■ Vp . (D.7) 

Also I calculate: 

^ ■ E ^Jkl ■ [Sku, ■ Fojl - ^liM • Kk\ = X ^jkl • ^kiM • Kl 
^ k,l=l k,l=l 



^Note that: J2k=i ' '^nmk = 6jn • 6pm - 6jm ' Spn for j,p, m,n = 1, 2, 3. 
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3 3 3 3 

P P P P P P P 

k,l=l m=l k,l=l m,n=l 

k,l=l ^in,n,r=l 

3 / 3 

= X! ^i'fc ■ ( X! ^nmfc ■ Sn " ^mZr ' -Bp^ — Qk ' E 
k,l=l m,n,r=l 

3 

= X! ['^jn ■ ^Im — ^jm " ^li^ ' Sji ' £mlr ' -Bp,, + {qAEp)j 
l,m,n,r=l 
3 

= -Yl ejnr ■ Sn " 5p,. + A Ep), = {Bp A s), + (g A Ep), , (j = 1, 2, 3) (D.8) 

n,r=l 

from which follows 

■ ^1 - ■ f1 ^sABp + {Vp A A Ep , (D.9) 

where the expression on the rhs of the arrow denotes the corresponding spatial part. With 
(D.7),(D.9) the spatial part of the rhs of (5.5c) reads as: 



— ■ i^^-lp -s^ -Vp 
m 2 ^ 



{Ep+VpABp)+vl-Ep-s 



^ ■ll-s^-Vp-v^^-Ep-Vp + ^-- A fip) - f ■ . Ep ■ Vp] , (D.IO) 



2 ip p p p p • 2 

and the spatial part of the Ihs of (5.5c) reads as: 

Ip-s'. (D.ll) 

Thus I have shown that the spatial part of equation (5.5c) is equivalent to equation (5.12c). 
Using (3.6) one also finds that the temporal part of equation (5.5c) follows from equation (5.12c) 
so that equation (5.5c) is equivalent to (5.12c). 



Appendix E 

In this Appendix I show that (6.8) is equivalent to (5.5). From (6. 2), (6. 7) follows 



H P ~ P P ^„P ~ J P ^„P P '~ P P P 

-Sux ■ S^p ■ 6^a] -Tp -U^ ■ F^^ = 5up -U^ -Tp ■ F^^ - -T^ - Up + -U^ -T^ , 

P P P P ^ P ~P PP~P PP 

^UOL ' ^ap ~ ^yp ' ' ' ^olI3 ~ ^pa ' ' '^v ~^ ^poi ' ' i 

P P P P P P ~ P ~ P P 

sL-F:, = 2.U:-T; -F;;^, {f,,u,p = l,...,A) (E.l) 

Note that: Eai^Pj ■ £apXa = <5i>p • Sf}\ ■ dja + Sua ■ 5f3p ■ + 5u\ ■ Spa- ■ Sjp 

■ Sup ■ Spa ■ S.yx - Sua ■ Spx ■ S^p - 5u\ ■ Spp ■ Sya for ly, (3, 7, p,\,a ^ 1, 4. 
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so that one gets by (5. 5c), (6. 1-2): 





i 

~ ~2 ' 




i 






i 








-f" 

va 


P P P 


i 

~2 




^ 4-m ^ 



2 

ap ' u! va ' OLp ' ui ' ^ ' '^p ' Lua ' a J 

p p p ~P P P P ~P P P P 

■ I ■ J J ■ J J A- P ■ J J ■ T ■ JJ A- P ■ I ■ J J ■ J J 

p p p p p p p p 

^pa/S-y ■ -T^ ■ ' + EpafS^ ■ F^^ ■ ■ ■ 

+£vaf3'y ■ Fp^ ■ ■ Up ■ - Eua/B-y " F^^ ■ ■ Tp ■ U^) — ■ Eypp^ ■ ■ ■ ■ 



■ Epupuj ■ Epaj3-f ■ [ Fp^ ■ -Uj^ 'U^ + Fp^ ■ -T^ 'U^\ 



Q ■ m 

^■9 ip'' t'' jf Tf - p"" jf t'' Tf] 

8 ■ m ' ''^ ' '^"l^'y ' I- /37 " ■ "^p ' ^uj ^13"/ ' '^a ' -'- p ■ "-^ajj 

6 rj,P jjP pP jjP 
Zi ' lib 

e • 9 IJ^P rj.P jjP jf ,-pP jf rj.P jjP^ 

4 m ' '"^^^ ' P"!^"^ ' I- P'y a ' V ' u} ^ ^ fS-y ' '^a ' -'^ v ' J 

6 P P P 



2 ■ m 

_ ^-9 ^ ^ jf jf ^_ p p t"" jf p'' jf 

~ A . m ' P'^t^'^ ' ^P<^P^ ' ^ f3j ' ^ a ' ^v ' '-^w O . rn ' ^"PP^ ' '^fpPl ' /3 ' ^7 ' ^ u)a ' ^ a 

6 ■ g 

= ^ — — ■ [{^va ■ ■ ^1^7 + ^f^"/ ' ' ^a;/3 + ^vf3 ' ^^J.■y ' ^uja ~ ^va ' ^fi-y ' ^uif3 ~ ^vy ' ^pf3 ' ^Loa 

p p p p G p p p p 

-5^/3 ■ Sf,a ■ ^u-i) ■ Fp^-U^ -T^ ' ^ ~ m ' ^^^'^ ' '^'^'^ " ' ^""^"^ '^'^ '^^ ' ^"""^ ' 

= ^ ■ K ■ u: ■ t; ■ < + f;„ . r:i + i- . f;„ ■ t; ■ u: ■ c/; 

Also one gets from (E.l): 

P P P P P P P P P P P ~ P P P P P P ~ P 

Sv.-dpF^^ + Up-S^^-U, .d,F^^ = 2-U^ -Tp ■ dpF^^ + 2 ■ U ^ ■ -T^-U, -d.F^^, 

P P P P P P P P P P ~ P P P P ~ P 

S,v-U^-U, .d,F,^ = Up -U^ -Tp -U, -d.F^^ + T, ■ U, -d^Fp,, 

p p p p p p p p ~ p p ~ p p 
p . Q ■ P ■ U = P ■ U ■ U ■ 1 ■ P 4- P ■ P ■ 1 

pv ^vlli u)p ^ p fiv ^v ^a (3 ajS ' fj.v va a ' 

p p p P P P P P ~ P 

Fpv ' Uy ■ S^p ■ Fp^ = 2 ■ Fp^ ' ■ ■ Tp ■ F^p , 

Suv ■ F^p ■ F ■ = F -U^ -U^ -Tn ■ F^a - F ■ ■ ■ ■ F^. . (/i = 1, •••,4) 



(E.3) 



Combining (5. 5), (6. 8), (E. 2-3) one observes that (6.8) is equivalent to (5.5). 
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Note that: Syp^,^ ■ e^pp^ = Spp ■ - Sp^ ■ S^p for ^j,,uj,(3,j = 1, ...,4. 
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Remarks accompanying the text 



[a]: For a particle with 'normal' intrinsic magnetic dipole moment one has by definition: g = 2. 



[b]: Transformations have been given in ||B162| , PS7CI|| and in ||CaM55|| for special cases and 
Jackson mentions the 'Pauli reduction' ||J ac76|| . The general transformation is straightforward 
and has been carried out by the author. It will be published separately ||Hei | . Note that 
it is consistent with the (semi-relativistic) Foldy-Wouthuysen transformation | )BD64| , |(JolV195| , 
|FW50| , [Fol62|| , because expanding the Hamiltonian H^^j^ up to second order in 1/c yields by 
(0.1) for g = 2: 



m 



-7? + - ■ 6 

M.op M.op/ 2 TM.op 



2 ■ m 



6 -|- ~* 6 ■ 

' ' ^op ' Bm,op ~^ ^ _ ^2 ' ^op ' ^^M,op ' ' ^M.op 



A E^j ^ ) + hermitian conjugate . 



The Darwin term does not appear here because of the vacuum Maxwell equations (see section 
2). Expanding Hj^^^^ only up to first order in 1/c yields for g = 2: 



H 



M.op 



m ■ f3 + 



2 ■ m 



(3-n 



■IT,, + e ■ (p., 

M,op M,op T^M.op 



m 



-B 



M,op ' 



which has the form of the Schroedinger-Pauli Hamiltonian. 



[c] : Thus the DK equations turn out to be Poincare covariant but not manifestly Poincare 
covariant whereas the Frenkel equations are manifestly covariant. 

[d] : The particle described in this work has arbitrary but nonvanishing charge, arbitrary in- 
trinsic magnetic dipole moment and vanishing intrinsic electric dipole moment. For a neutral 
particle the equations can be easily modified. 

[e] : Here Eijk is the antisymmetric symbol with £123 = 1 and S denotes the Kronecker delta. 
All three-component quantities r^^ , , a, . . . denoted by an arrow are column vectors. The 
components aj of any a are defined by: 

a = (ai, 02, a^y . 

The transpose of a three-component quantity is denoted by 'f. Therefore rl^^pl^^a^ . . . are 
row vectors. 

[/]: In this paper the nabla operator V always acts on functions depending on fj^j,t,p^j,a 
and it is the gradient w.r.t. f^^. 

[g]: Note that in this paper the multiplication symbol '■' always denotes matra multiplication 
and that a single number is a 1 x 1 matrix. 

To avoid the mushrooming of the bracket symbol I avoid its use even in places where it 
would usually help to find the correct order of matrix multiplications to be performed. If for 
example a matrix product like 

TT ^ ■ R ■ 7? 

"m -^m "m 
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occurs (in this example associativity of multiplication does not hold!), then the matrix struc- 
tures of the factors suggest the correct order of the matrix multiplications. In the present 
example one has: 

TT ^ ■ B ■ TT = (tT^ ■ B ) • TT 

because 

is a meaningless expression. Note also that the multiplication symbol 'A' always the denotes 
the vector product of three-component quantities. 

[h]: The nabla operator Vp in this paper always acts on functions depending on rp,t,Vp,s 
and it is the gradient w.r.t. fp. My notation is chosen so as to indicate that the functional 
dependence of Ep on rp,t is the same as the dependence of Ej^^ on rj^.j,t and likewise for 5^^. 
The explicit way in which the 'P' fields are derived from the 'M' fields is shown in (B.21). 



[i]: I am dealing with special relativistic space-time positions, tensors, pseudotensors, tensor 
fields and pseudotensor fields and I do so by using the complex notation where the fourth (=tem- 
poral) component is imaginary. I use this convention following the usage in much of the litera- 
ture on the relativistic SG force [^orS|, [Fri2^ , |Goo62| , |Nyb62| , |iNyb64 |Pla664 |Pla66b| , |Raf7D[ . 
Note also that Einstein's summation conventions are applied to Greek indices. The Greek in- 
dices assume the values 1,2,3,4. One of the advantages of the complex convention is that only 
lower components occur. For textbooks using this convention, see for example [ [Moe72| , |Syn58| . 
Of course one could use the real convention where, however, one has to distinguish between 
covariant and contravariant components. Most textbooks use the real convention, e.g. ||Jac75 



[j]: Any antisymmetric tensor of rank 2 can be characterized by two three-component 
quantities b resp. d, called the spatial resp. temporal part of the tensor, and they are defined 
by: 

1 3 

= 9 ■ H ^Jki ■ Nki , dj = . {j = 1, 2, 3) 

k,i=i 

I denote this correspondence as follows: 

N ^ (6, d) . 

Thus the spatial resp. temporal part of M is given by flp resp. i ■ Sp, i.e. 

M ^ iflp,i ■ Sp) . 

[k\. The constraint (3.8b) follows from (3. 2), (3. 6) and states that the particle has no intrinsic 
electric dipole moment ||Cor68| , [Fre26| , |Nyb64|| . 



[/]: In this paper the partial derivatives (9f,...,94 always act on functions depending on 

P P P P 
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[m] : One also gets this Hamiltonian if one neglects terms of second order in spin in the Hamil- 
tonian given by Corben. For details, see ||Cor68|| , especially Chapter 7 thereof. If one neglects 
all spin terms then H reduces to the following well known expression: 



H 



1 „p „p m 

n„ ■ n„ + — 

2-m ^ ^ 2 



See for example ||Ba64| , |Gol80| , |Jac75| 



[n]: The assumptions made on are as follows. The proper time dependence in the SG term 
on the rhs of (5.1b) is assumed to come in only via the proper time dependence of X ,U , S . 

P P 

To get a useful class of allowed Y , I assume in addition that the X dependence only comes in 

P P ~ 

via F and its space-time derivatives. Specifically I assume that Y^ 
following arguments: m,Ui , f/4 , 

P ■ 

the six components of F . This function is supposed to be a polynomial in all its arguments 
(except m) and first order in spin. Thus the coefficients of this polynomial are functions of m 
and it turns out that they are just powers of m times dimensionless numbers. 



^, Y[ are functions of the 

P P ^ 

Ua , the six components of S and all space time derivatives of 



[o\. Another remark on the quantum mechanical aspect is in order. In deriving the classical 
equations from the Dirac equation (plus the Pauli term) one does not get a unique answer 
because one depends on the choice of the proper operators. For example the Frenkel equations 
can be obtained by using a certain special relativistic generalization of the Foldy-Wouthuysen 
transformation ||B162| , pS70| , pS72| , [Hei]| with the emphasis on the Newton- Wigner position op- 
erator whereas the GNR equations can be derived by the 'Gordon decomposition' |[Raf 70|| . Note 
again that in the present work I am only interested in the classical aspect. Further details on 
the quantum mechanics including the opinions of Pauli and Bohr are discussed for example in 
pSV86|,P^,^eT72|. 



[p\. The Lorentz group consists of the homogeneous part of the Poincare group, i.e. it does not 
contain the space-time translations. Note that the Poincare group is also called the inhomoge- 
neous Lorentz group. A Poincare transformation is composed of a translation a and a Lorentz 

p 

transformation L so that the space-time position X transforms as: 



(/i=l,...,4) 



Since L belongs to the Lorentz group it satisfies 



L 



6. 



up 



1, 



.4) 



The restricted Poincare group (=proper orthochronous Poincare group) contains those Poincare 
transformations, where: 



det(L) 



L44 > . 



Hence the restricted Poincare group contains neither the parity transformation nor the time 
reversal transformation. For more details on the subgroups of the Lorentz group resp. Poincare 
group, see for example PLT75| , pW89|| . 



[q\. This follows because the relations (3.1), (3. 6), (3. 9) between the variables X^,S'',F^ and 
the variables , t, s, Bp , Ep are the same for every Poincare frame. 
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[r] : If one performs a Poincare transformation, then (as shown in subsection 5.3) the equations of 
motion (2.11) for the 'P' variables remain the same (except that the 'P' fields have transformed 
in a specified way). Under the same Poincare transformation the equations of motion (1.5) for 
the 'M' variables also transform in a definite way. The transformed equations of motion for the 
'M' variables can be derived from the transformed equations of motion for the 'P' variables in 
the same way as the original equations of motion (1.5) for the 'M' variables were derived in 
section 2 from the original equations of motion (2.11) for the 'P' variables because the relations 
(2.1), (2. 3a), (2. 4-5), (2. 7) between the 'M' variables and the 'P' variables are the same for every 
Poincare frame. The result is that the transformed equations of motion for the M' variables 
are the same as the original equations of motion (1.5) for the 'M' variables (except that the 
'M' fields have transformed in a specified way). Therefore the DK equations (1.5) are Poincare 
covariant. 
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